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ABSTRACT 

In this paper, the second of a series, we study the stellar dynamical and evolutionary processes 
leading to the formation of compact binaries containing neutron stars (NSs) in dense globular 
clusters. For this study, 70 dense clusters were simulated independently, with a total stellar 
mass ~2x 10 7 M@, exceeding the total mass of all dense globular clusters in our Galaxy. 

We find that, in order to reproduce the empirically derived formation rate of low-mass 
X-ray binaries (LMXBs), we must assume that NSs can be formed via electron-capture super- 
novae with typical natal kicks smaller than in core-collapse supernovae. Our results explain 
the observed dependence of the number of LMXBs on "collision number" as well as the large 
scatter observed between different globular clusters. We predict that the number of quiescent 
LMXBs in different clusters should not have a strong metallicity dependence. 

We compare the results obtained from our simulations with the observed population of 
millisecond pulsars (MSPs). We find that in our cluster model the following mass-gaining 
events create populations of MSPs that do not match the observations (either they are incon- 
sistent with the observed LMXB production rates, or the inferred binary periods or companion 
masses are not observed among radio bMSPs): (i) accretion during a common envelope event 
with a NS formed through electron capture supernovae, and (ii) mass transfer (MT) from a 
WD donor. Some processes lead only to a mild recycling - physical collisions or MT in a 
post-accretion induced collapse system. In addition, for MSPs, we distinguish low-magnetic- 
field (long-lived) and high-magnetic -field (short-lived) populations, where in the latter NSs 
are formed as a result of accretion induced collapse or merger induced collapse. With this 
distinction and by considering only those mass-gaining events that appear to lead to NS re- 
cycling, we obtain good agreement of our models with the numbers and characteristics of 
observed MSPs in 47 Tuc and Terzan 5, as well as with the cumulative statistics for MSPs 
detected in globular clusters of different dynamical properties. 

We find that significant production of merging double NSs potentially detectable as short 
7-ray bursts occurs only in very dense, most likely core-collapsed clusters. 

Key words: stellar dynamics - binaries: close - binaries: general - stars: neutron - pulsars: 
general - globular clusters: general - X-rays: binaries. 



1 INTRODUCTION 

Neutron stars (NSs) are seen in globular clusters (GCs) via their 
(current or past) participation in interacting binary systems. From 
the earliest observations of X-ray binaries in GCs it has been noted 
that their abundance per unit mass is ~100 times greater in GCs 
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than in the Galaxy as a whole dKatzlll975t) . This was understood 
to be a consequence of the high stellar density of GCs, which may 
lead to the creati on of compact NS binaries in close stellar encoun- 
ters dClark|[l975l) . NS binaries are often relatively easy to identify 
in GCs (compared to other binaries), and they allow us to constrain 
the dynamical history of clusters and the evolution of NS properties 
in binary systems. They have been observed so far in three guises: 
(i) as bright low-mass X-ray binaries (LMXBs); (ii) as quiescent 
LMXBs (qLMXBs); and (iii) as binary or single millisecond pul- 
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sars (MSPs). In the future, they could also be detected as sources 
of gravitational radiation. 

Bright LMXBs have X-ray luminosities of typically 10 35 - 
10 38 ergs s _1 . In bright LMXBs, the NS is actively accreting 
from a companion star via Roche-lobe overflow and an accretion 
disk. Thirteen sources are known in Galactic GCs, of which 12 
are known to contain NSs through the detection of thermonuclear 
bursts. Six are transient systems (two have remained in outburst 
for more than a decade; the others have outburst time-scales of 
weeks), while the other seven seem to be persis tent. Orbital peri- 
ods have now been identified for seven systems (jVerbunt & Lewinl 
120061 : iDieball et al]|2005l ; lAltamirano et alj|2007l) . of which three 
are below one hour, indicating ultracompact systems composed of 
NSs with white dwarf co mpanions are common in globular clus- 
ters dDeutsch et al.ll200oh . Efficient surveys of the X-ray sky give 
us confidence that we have detected all bright LMXBs that have 
existed in Galactic GCs during the past 10 years. 

Quiescent LMXBs have lower X-ray luminosities (10 31 < 
Lx < 10 34 ergs/s), indicating that accretion is not active or is 
significantly reduced during quiescence. Their X-ray luminosity 
is thought to derive at least in part from the re-radiation of heat 
from t he core of the NS, w hich was generated during accretion out- 
bursts terown et alj|l998[). though other mechanisms are also sug- 
gested (see, e.g., ICampana et al.l d 19981) ). They can be identified 



through their positional coincidence with the location of a bright 
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transient LMXB in outburst (Wijnands et al. 2005). and/or through 

their characteristically soft X-ray spectrum iRutledge et alj|2002l) . 
The total number of qLMXBs in the Ga lactic GC system is in- 
ferred to be ~ 100-200 drleinke et ai1l2003l) . although only a frac- 
tion of the Galactic clusters have been studied sufficiently to iden- 
tify qLMXBs. 

As of this writing 129 millisecond radio pulsars have been ob- 
serve d in GCs, with mo re than 1 0, 20 and 30 seen fro m the clusters 
M28 dStairs et al.l2006l) . 47 Tuc dCamilo et alj2000l) . and Terzan 5 
dRansom et al.l 120051) , respectively These millisecond pulsars 
(MSPs) are thought to be descendants of LMXBs (for reviews se e 
iBhattacharva & van den Heuvel|[l99~il IPhinnev & Kulkarni|[l994T) . 
It is estimated t hat of order 1000 M SPs are present in the Galac- 
tic GC system drleinke et all l2005bl) . of which half or more are 
potentially detectable. Current pulsar searches have reached the 
bottom of the pulsar luminosity function in very few clusters 
(e.g., 47 Tuc and M15). As an example, the total radio luminos- 
ity f rom Terzan 5 suggested that this cluster may contain 100 pul- 
sars dFruchter & G oss 1995) — this number is 3 times more than the 
number of identified pulsars. 

Since NSs in GCs appear to be excellent tests for our un- 
derstanding of both stellar dynamics and neutron star evolution, 
it is not surprising that this issue has attracted significant contin- 
uing attention. The largest challenge for these studies is computa- 
tional. It has been estimated that a cluster consisting of 10 6 stars 
(~ 250, OOOMq) will contain only a few hundred NSs, and mod- 
els with a smaller number of stars will not therefore have suffi- 
cient statistics. Direct TV-body met hods are un likely to address the 
million-body problem before 2020 JHu3|2006|) . Monte Carlo meth- 
ods now allow million-bo dy simulation s for spherically symmetric 
systems with single stars dGiersj|2006h and with large fractions of 
primordial binaries, including direct integrations of all binary en- 



1 Not all qLMXBs show the characteristic soft X-ray spectrum 
IWiinands et alj2005h . though it seems that the majority do. 

2 See http://www.naic. edu~pfreire/GCpsr.html for an updated list. 



counters dFregeau & Rasioll2007l) . but a full impleme ntation of bi- 
nary star evolution is still lacking dRasio et"al]|2006l) . As a result, 
most previous theoretical stu dies of NSs in clu sters were performed 
using analytic esimates (e.g. jHut et aljfl99ll) . or no binaries (e.g., 
IPrukieJ 19961) , or simple scattering e xperiments, with a static back- 
ground of single or binary stars (e.g..lSigurdsson & Phinnevll 19951 : 
iRasio et al.ll2000l : iKuranov & Postnovll2006l) . We chose instead to 
use a method that treats in detail the stellar evolution of a large 
population of stars and binaries in a dense cluster environment, but 
with a highly simplified description of the stellar dynamics which 
is not fully self-consistent. 

In this paper we present the results of our ongoing study of 
compact binary formation in GCs, concentrating on the forma- 
tion of binaries with NSs. For the first time we have used cluster 
simulations that include a dynamically evolving system and rep- 
resent consistently an entire stellar population with numbers of 
stars and binaries comparable to those in observed GCs. In addi- 
tion, several statistically independent simulations were performed 
for each cluster model, to show how different the outcomes can be 
with slightly different but equivalent initial conditions. Some pre- 
liminary results o n NS binaries from our work were reported in 
llvanova & Rasiol d2004l); llyanova et al.l d2005bl) ; llvanova & Rasiol 
d2005d) ;l Rasio et al .1 d2006l) . and detailed results for the formation 
of bin aries containing whi te dwarfs (WDs) were presented in Pa- 
per I dlvanova et al.l 20061) . We refer the reader to Paper I for a 
complete description of the physical processes of formation and 
destruction of mass-transferring compact binaries, as well as our 
numerical methods. We outline updates to the method and our as- 
sumptions that are especially relevant to NSs in Section 2. Section 
3 outlines the scenarios for formation of MSPs in GCs. In Section 4 
we describe the range of GC models that we consider. Formation 
and retention of NSs are analyzed in Section 5, and the follow- 
ing section is devoted to the formation of binaries containing NSs. 
We discuss the formation of LMXBs in Section 7, and MSPs are 
considered in Section 8. Section 9 briefly describes double neutron 
stars. We conclude by summarizing the connection between our re- 
sults and the observations, as well as giving predictions for future 
observations and discussing constraints on NS formation. 



2 METHOD 

Our numerica l methods gen erally remain the same as described in 
llvanova et alj |2005a, 2006). To study the population of NSs, we 
introduced a number of updates to the population synthesis model 
specifically for NS formation, and we revised the treatment of dy- 
namical events with NSs. We also optimized numerical aspects of 
our code for more efficient use of parallel supercomputers. 



2.1 Population synthesis updates 

We adopt the binary evol ution model from the populat ion syn- 
thesis code StarTrack dBelczvnski et alj I2002L 120081). In our 
first s tudy of binary populations in GC cores dlvanova et alj 
l2005al) we used a previo us version of StarTrack described 



_ Belczvnski et alj d2002l) In ou r study of close binaries with 
white dwarfs ( llvanova et ai] |2006l , paper I) we employed a more 
recent version of St arTrack available at the time (as in 
iBelczvnski et alj|2008l . but including updates only for WD evo- 
lution compared to the previous version). In this paper we in- 
corporated all th e latest updates of StarTrack, as described in 
IBelczvnski et alj d2008l) . Below we outline several of the most 
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important changes in StarTrack that affect the formation and 
evolut ion of NSs compared t o our previou s studie s of N Ss in 
GCs dlvanova & RasTol |2004 llvanova et al.1 l2005bl 120061) . We 
also provide detailed descriptions of fu rther modification s to th e 
StarTrack model, beyond those in iBelczvnski et alj J2008h . 
which we developed for the treatments of electron-capture super- 
novae and common-envelope events. 



2.1.1 Natal kicks and supernova disruptions 

At the moment of formation, both NSs and BHs receive addi- 
tional speed (a natal kick), most likely due to asymmetric su- 
pernova ejecta dFrverj|2004 . Although most recent simulations 
are in relatively good agreement with the measured distribution 
of pulsar velocities , the agreement is not yet firmly established 
dScheck et al.ll2006j) . In this paper we adopt t he most recently de - 
rived pulsar kick velocity distribution from iHobbs et al.l d2005h . 
which is a Maxwellian distribution with one-dimensional RMS 
velocity o — 265 km s _1 (the mean three-dimensional velocity 
is ~ 400 km s~ ). In c ontrast to some earlier studies (e.g., see 
lArzoumanian et al]|2002h . no evidence for a bimodal velocity dis- 
tribution is present. 

At the moment of the explosion, the supernova progenitor can 
be a binary companion. Because of the kick, the binary can become 
unbound (and the binary components will separate) or the entire 
binary can be affected by the natal kick. In our previous studies, in 
order to find the velocities of the component s of an unbound binary, 
we used the derivation bv lTauris & Takerisl dl998h . where the pre- 
supernova orbit is assumed to be circular. The approach used in 
the current version of StarTrack allows us to properly calculate 
cases with an arbi trary pre-explosion orb ital eccentricity (for full 
details, see §6.3 in lBelczvnski et al.l2008l) . 



2.1.2 Electron-capture supernovae 

Through our studies of accreting WDs dlvanova et alj l2006r) we 
found that electron-capture supernovae (ECSe) could be the dom- 
inant source of reta ined NSs in clusters. Simil ar results were also 
recently obtained in lKuranov & Postnovl d2006h . We have therefore 
considered in great detail all possible types of ECSe for this work. 

It has been argued that when a degenerate ONeMg core 
reaches M ccs = 1.38Mq, its collapse is triggered by electron cap- 
ture on 24 Mg and 20 Ne before neon and subse quent burnings start 
and, therefore, before t he iron core formation dMivaii et al.ll 19801; 
Nomotd 1 1984 1 19871 : iTimmes & Wooslevl 1 19921 ; ITimmes et al.1 



1994 . The explosion energy of such an event i s significantly lower 



than that inferred fo r core-collapse supernovae dDessart et al. 120061 : 
iKitaura et ai]|2006h . There are several possible situations when a 
star can develop a degenerate ONeMg core that will eventually 
reach M ocs . 

First, this can occur during the normal evolutio n of single 
stars. It was initially proposed by iBarkat et al.1 jl974l) that a star 
of 7-10 Mq, after non-explosive carbon burning, develops an 
ONeMg core. If the initial core mass is less than required for 
the neon ignition, 1.37M0, the core becomes strongly degener- 
ate. Through the continuing He shell burning, this core grows to 
M ocs . In more massive stars, > 10 Mq, carbon, oxygen, neon and 
silicon burnings progress under non-degenerate conditions, and, in 
less massive stars, ONeMg cores never form. The formation of de- 
generate or non-degenerate ONeMg cores depends on the He core 
mass at the start of the asymptotic giant branch, where ~ 2.5 Mq 



is a rough boundary between the cases dNomotd[l984l) . If the ini- 
tial He core mass is below 1.83 Mq, no off-centre ignition will 
happen — the carbon core burning will occur when the degenerate 
core reach es the Chandrasekh ar mass, resulting in a thermonuclear 
explosion dHurlev et alfeOOd. also J. Eldridge 2006, priv. comm.). 

The conditions for E CS were shown to occur in single stars of 
initial mass 8-10 Mq by |Nomotd l fl984h . This critical mass range 
depends on the properties of the He and CO cores, which, in turn, 
are highly dependent on the mixing prescription (semiconvection, 
overshooting, rotational mixing, etc.) as well as on the adopted 
opacities. As a result, the range varies between different evolu- 
tionar y codes (see discussion in IPodsiadlowski et alj|2004 ISiessI 
l2006h. and i s generally re duced compared to that proposed initially 
bv lNomotddl984ll987l) . In the code that we use for our cluster 
simulations, a non-degenerate ONeM g core is formed when the ini- 
tial H e core mass is about 2.25 Mq dPols et al.lll998l ; lHurlev et al.1 
l200d) and the range of initial masses for single stars of solar metal- 
licity that leads to the formation of such a core is 7.66 to 8.26 
Mq. The equivalent mass ranges are from 6.85 to 7.57 Mq and 
from 6.17 to 6.76 Mq for single stars with GC-like metallicities, 
Z = 0.005 and Z = 0.0005, respectively. 

The initial stellar mass is not the only parameter that defines 
whether a star will form a degenerate ON eMg core. It was recently 
pointed out bvlPodsiadlows ki et al.ld2004h that the range of progen- 
itor masses for which an ECS can occur depends also on the mass 
transfer history of the star, and therefore may be different for binary 
stars, making it possible for more massive progenitors to collapse 
via ECS. We will refer to the single and binary scenarios for ECS 
in non-degenerate stars, as described above, as evolution-induced 
collapse (EIC). 

The second possibility for an ECS to occur is by accretion on 
to a degenerate ONeMg WD in a binary: accretion-induced col- 
lapse (AIC). In this case, a massive ONeMg WD steadily accumu- 
lates mass until it reaches the critical mass M cca ( for more detail on 
the adopted model of the accr etion on the WD see llvanova & Taaml 
|2004 IBelczvnski et alj|2008l) . In addition, ONeMg WDs can be 
formed from CO WDs via off-centre carbon ignition if the mass 
of CO WDs is above 1.07 Mq and the accr etion rate is high, 
M > 2.7 x 10~ 6 Moyr _1 dKawai et al]|l987t) so that eventually 
this also leads to AIC. 

The third case that we consider is coalescing double WDs with 
a total mass exceeding M ecs . To distinguish it from AIC, we will 
refer to this case as merger-induced collapse (MIC). The nature of 
NS formation is the same as in the case of AIC (accumulation of 
mass by a massive ONeMG WD). In the case of coalescing CO 
WDs, then, as in the case of single stars, off-centre carbon ignition 
occurs first, which quiescently conv erts the star into an ONeM g 
core, and this proceeds with an ECS dSaio &Nomotdll985ll2004l) . 
Although the approach above was proposed for coalescing WDs 
in binaries only, we assume that in the case of collisions between 
two WDs it is applicable as well: except for the rare cases of near 
head-on collisions, a less massive WD will be tidally disrupted and 
accreted at high rate on to the more massive WD. 

We therefore assume in our simulations that a NS will be 
formed via ECS in the following cases: 

• if at the start of the AGB the He core mass is 

1.83 Mq < Mc.BAGB < 2.25 Mq (EIC) 

• if an accreting ONeMg WD reaches M 8C s (AIC) 

• if the total mass of two merging or colliding WDs exceeds 
Mecs (MIC). 

For NSs formed via ECS we assume that the accompanying 
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natal kick is 10 times smaller than in the CC case. This assump- 
tion follows from the results of numerical simulations which find 
that the SASI instability, required by current understanding for the 
large explosion asymmetry in the case of core-collapse supernovae, 
fails to develop and resulting kick s are significantly smaller and on 
overa ll do not exceed 100 km/s l lBuras et"al]|2006l . iKitaura et al .1 
120061) . We will provide the separate statistics for NSs created by 
different channels (core-collapse, EIC, AIC or MIC), so that if one 
of the channels later turns out to have been overestimated it will be 
easy to recalibrate. 

The gravitational mass of the newly formed NS is 1.26 Mq. 
This is ~ 0.9 of its baryonic mass (which cannot exceed the pre- 
collapse cor e mass) after subtract ing the binding energy and is cal- 
culated as in lTimmes et al.l l ll996l) : 



Afbaryon _ Mg 



AM = 0.075AC 



(1) 



where the masses are in solar masses. As the angular momentum of 
the binary is conserved, in the case of no kick, the binary widens. 

The condition for ECS to occur depends essentially on the 
central density of the object. As a consequence, a rapidly rotating 
WD can reach a much higher mass than the Chandrasekhar limit 
before the central density becomes high enough for electron cap- 
tures on 24 Mg and 20 Ne to occur. For example, its mass can be 
as high as 1.92 Mq if the ratio of the WD rotational energy to 
the gra vitational binding e nergy is 0.0833 and the star is highly 
oblate l lDessart et alj|2003) . Such rapidly rotating heavy WDs can 
be formed, e.g., during the coalescence of two WDs with total mass 
exceeding the Chandrasekhar mass. During such a merger, less than 
0.5 p er cent of mass will be lost from the system dGuerrero et al.l 
120041) . The collapse of a rapidly rotating WD can therefore lead to 
the formation of a more massive and very fast spinning NS. 



2.1.3 Convective and radiative envelopes 

The presence or absence of a deep outer convective envelope influ- 
ences a star's behavior during various phases of close binary inter- 
actions (e.g., angular momentum loss via magnetic braking, tidal 
interactions, stable or unstable ma ss transfer). In particular, it has 
been pointed out in llvanov3 d2006h that the absence of a deep outer 
convective zone in low-metallicity GCs may expl ain the preferen- 
tial formation of LMXBs in metal-ric h clusters dGrindlavlll993l ; 
iBellazzini et aljl995tlzepf et alj2006l) . 

In this study, for main-sequence stars we adopt the metallicity- 
dependent ma ss range to develop a c onvective envelope described 
by (eq.(10) in lBelczvnski et alj [2008) In addition, in our previous 
work, all giant-like stars were assumed to have outer convective 
envelopes. Now stars crossing the Hertzsprung gap and stars in the 
blue loop are not assumed to have convective envelopes if their 
effective temperature satisfies log 10 T c g > 3.73. 



2.1.4 Common-envelope phases 

For com mon-envelope (CE) ev ents we use the standard energy for- 
malism l lLivio & Sokerl fl988l) . In this approach, the outcome of 
the CE phase depends on the adopted efficiency for orbital energy 
transfer into envelope expansion energy a cc and on the donor en- 
velope central concentration parameter A. The commonly used pre- 
scription is a cc x A = 1, and this is adopted in our models as well. 
However, for He stars with an outer convective envelope we deter- 
mi ne A from the set of d etailed He star evolution models calculated 
bv llvanova et~alld2003l) . In this case. A = 0.3i?~ ' 8 , where R is the 



radius of the He star in solar radii. As many He stars end their lives 
as NSs, the CE evolution of He stars is an important channel for 
the formation of close binaries with NSs. Indeed, the prescription 
above was recently used in a stu dy of double neutron sta r formation 
as short 7-ray burst progenitors (Belczvnski et al. 20o3). 



2.2 Dynamical events 

2.2.1 Binary formation via physical collisions with giants 

We consider in our models the possiblity of (eccentric) binary for- 
mation via physical collisions involving a red giant (RG). For bina- 
ries formed through NS-RG collisions, the final binary separation 
at and eccentricity et depend on the closest approach distance p 
dLombardi e7ai]|2006l) and can be estimated using the results of 
hydrodynamic calculations as 



e f = 0.88 



3Rrg 



Of 



3.3(1 



(2) 



(3) 



More details on our treatment of these collisions can be found in 
llvanova et alj J2006h . In this paper, we refer to this event as to a 
dynamical common envelope (DCE). 



2.2.2 Collisions and mergers with NSs 

We assume that the amount of matter accreted by a NS following a 
merger or collision depends on the evolutionary stage of the other 
star. If a NS collides or merges with a MS o r a He MS star, we 
assume that it accretes no more than 0.2 Mq dRosswo 120061) . In 
the case of a collis ion with a giant-like st ar, the maximum accreted 
mass is 0.01 Mq dLombardi etai]|2006l) . When the other star is a 
WD, all the WD mass is accreted. In essentially all mergers and 
collisions, the accreted mass must have gone through an accretion 
disk, implying that, if a significant amount of matter is accreted, a 
recycled pulsar will be formed. 



2.2.3 Triples and their treatment 

Stable hierarchical triples can be formed through binary-binary en- 
counters. Although these triples would be stable in isolation, it is 
likely that they will be destroyed during their next dynamical en- 
counter. As there are no developed population synthesis methods 
for handling triple star evolution, we cannot keep these triples in 
our simulations and instead we have to break them immediately af- 
ter formation into a binary and a single star. We do this while con- 
serving energy: the energy required to eject the outer companion is 
balanced by shrinking the inner binary orbit. The outer companion 
is released unless the required energy is such that the inner binary 
merges. In the latter case the inner system is allowed to merge and 
the outer companion is kept in a new, wider orbit to form the final 
binary system. 

It is possible that in a triple the inn er orbit's ecc entricity will 
be increased via the Kozai mechanism dKozai|[T962l) . This secu- 
lar coupling causes large variations in the eccentricity and incli- 
nation of the orbits and could drive the inner binary of the triple 
system to merge before the next dynamical interaction. The maxi- 
mum eccentricity in systems with large initial inclinations ip, such 
that sin ^0 > (2/5) 1/2 (ip > 39°) is (e.g., Ilnnanen et al lll997l: 
lEggleton & Kiseleva-Eggletonll200ll) : 
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e max ~ y/l-5/3 cos 2 (i ) . (4) 

The period of the cycle to achieve e max is dlnnanen et al Jl 19971 ; 
iMiller & Hamiltorj|2002h 

0.421n(l/ ei ) ( mi + m 2 b 3 , V /2 / feg V /2 
v/sin 2 (i ) -0.4 V m ° a ?/ V Gm °/ 

where ej is the initial eccentricity of the inner binary, mi, mi and 
m are the masses of the inner binary companions and the mass of 
the outer star, aj and a are initial orbital separations for the inner 
and outer orbits, and b = a (l — e ) 3//2 is the semiminor axis of 
the outer orbit. 

We compare the Kozai time-scale tk oz with the collision time 
Tcoii (computed as in llvanova et al . 20052): 



r coll = 8.5 x 10 12 yr P~ A/3 M^n^v^ 1 x 



(6) 



1+913 



(Aftot + (M)) 



Here P t( j is the triple period in days, M to t is the total triple mass 
in Mq, (M) is the mass of an average single star in M©, Wio = 
«oo/(10 km/s) and ns = n/(10 5 pc -3 ), where n is the stellar 
number density. 

If tkoz > T co ii, the Kozai mechanism does not affect the 
triple evolution before the next encounter occurs, and we break 
the triple as describe above. However, for triples with tkoz < 
Tcoii, we consider the circularization time-scale r c i rc (see §3.3 in 
iBelczvnski et alj|2008l) . The characteristic time-scale for the inner 
binary separation to decre ase through tidal dissipation is r c i rc /2e 2 
dMazeh & Shaharrtfl97^) . If r ciTC /2ef < r co ii, the inner binary 
will shrink until one of the components overfills its Roche lobe 
or its separation will sm all enough to cause the Darwin insta bil- 
ity resulting in a merger fcggleton & Kiseleva-Eg gleton 200lt). In 
the opposite case, when r c i rc /2e 2 > T co n, we assume that mass 
transfer starts if e max is such that at least one of the stars in the 
inner binary overfills its Roche lobe at pericentre. As a result, triple 
formation may enhance the formation of mass-transferring binaries 
with a NS. 

For sufficiently compact inner binaries, relativistic effects may 
play an important role. The relativistic precession period P m of the 
inner binary, to first post-Newtonian order, is (see I Weinberdl 1 9721 
p. 197): 



27rc 2 a 5/2 (l 



(7) 



3G 3 / 2 (mi +m 2 ) 3 / 2 

If tko z > P P r, Kozai cycles are strongly suppressed dHolman et al.l 
1 19971) . 

We expect that some dynamically formed triples can also un- 
dergo significant secular eccentricity evolution e ven for orbital in - 
clinations smaller than the Kozai ang le (see e.g. lFord et aT]|2000h . 
but here we neglect this possibility. 



3 RADIO PULSAR MODELING 
3.1 Slow, non-recycled radio pulsars 

The time over which a pulsar will spin down from its initial period 
Po to a current period P is 

(8) 



where PP = B 2 /10 39 [s]. With the assumption that the magnetic 
field B does not decay, PP remains constant with time. 

Several observed radio pulsars in GCs are slowly spinning 
(P > 0.1s) and have high Ps (higher than can be produced by 
gravitational acceleration in the cluster potential) implying rela- 
tively high B fields > 1 11 Gauss. The se include NGC 6 342A 
jvan Kerkwiik etalJl2000l), NG C 6624B jBiggs et al.l [l994h . and 
NGC 6440 A dLvne et alJI 19961) . Their inferred characteristic ages 
range from 10 7 yr to 2 x 10 s yr, compared to 10 11 — 10 13 yr for 
"normal" MSPs. These pulsars cannot have been formed in their 
current state through core-collapse SNe, which only occur within 
the first ten million years. The inferred birthrate for this class of 
objects is similar to that of "normal" MSPs, but observationally 
they seem to be concentrated in the densest clusters, and 2 of 3 are 
isolated pulsars. 

A strong magnetic field of > 10 11 -gauss can be produced 
by the collapse, assuming flux conservation, of a typical magnetic 
WD with B ~ 10 6 -gauss. Observed magnetic WDs typically have 
highe r masses than nonmagnetic WDs dWickramasinghe & Ferrariol 
120051) . indicating that they may dominate the population of WDs 
that undergo AIC. Stellar dynamo theory sup ports the link between 
strong magnetic fields and massive WDs ^Thompson & Duncan! 
1 1993b . It is also possible that a strong magnetic field may be gen- 
erated during the merger of a WD pair, where one or more WDs 
had a strong magnetic field, leading to the production of 10 12 -gauss 
magne tic fields, or even magnetar B field strengths of 10 14 -f- 10 15 - 
gauss jKing et allbOQll ; iLevan et alj|2006t IChapman et"aT1l2006h . 
MSPs formed with such high magnetic fields will have relatively 
short lifetimes, and thus will make up only a small portion of the 
observed globular cluster MSP population. In our simulations, we 
do not assign a specific magnetic field to formed neutron stars. 
However we keep a record of their formation and therefore can 
test whether the population of AIC and MIC NSs is formed at the 
required rates to explain the population of slow GC radio pulsars, 
assuming that their life-time is only 10 8 years. 

We refer to young, recently formed pulsars that have not yet 
been recycled as "la" and "lb" in Fig.Q] As a standard picture, we 
consider that the origin of a NS does not affect its further evolu- 
tion if a MT event occurred in a dynamically formed binary (c.f. 
post-AIC evolution), but we note that it is possible that MIC NSs 
might have a higher than usual magnetic field and their subsequent 
evolution might be different - they can be seen only for a shorter 
time after their formation as young slow pulsars and they might be 
not able to accrete and be spun up. 



3.2 Mild recycling in post-AIC systems 

Very little is known about whether a magnetic field is always re- 
duced during mass accretion, independently of its initial strength. 
E.g., the UCXB and 7.7 second X-ray p ulsar 4U 1626-67 ha s 
a strong (4 x 10 12 -gauss) magnetic field dOrlandini et al.|[l998h . 
On the other hand, known magnetic fields in accreting MSPs 
are usually 10 7 - 10 9 gauss (e . g., 1 - 5 x 10 s -gauss in XTE 



PP 



J1751-305, 15i Salvo & Burderil 420031) ; < 3 x 10 8 -gauss in IGR 
J00291+5934lTorres et aljj2007l) : < 10 9 -gaussinXTE J0929-314, 
iGallowav et aljfc002h ) and are consistent with the magnetic fields 
of MSPs. Even though 4U 1626-67 is evolving through the accre- 
tion phase and has passed the stage when the MT rate is maxi- 
mum, the NS has kept a high magnetic field and has not been 
spun up to millisecond periods. To explain this, it has been sug- 
ested that the NS in 4U 1626-67 was recently formed via AIC 
Yungelson et alj2002l) . Several radio pulsars in the Galaxy possess 
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relatively high B fields and lo ng periods, with lo w-mass compan- 
ions in low-eccentricity orbits dBreton et alj2007 t). They argue that 
these systems can be well-explained through AIC. 

It is unclear why post-AIC binaries may yet fail to produce 
MSPs. Reasons could include the high initial B field and mod- 
erately fast spin of newly-formed AIC NSs, which could disrupt 
spin-up; or simply that after sufficient accretion onto the WD to 
reach the Chandrasekhar limit, there is usually not enough remain- 
ing mass in the companion to spin the NS up to millisecond periods. 
We will show below that producing MSPs in the same binaries in 
which AIC occurred contradicts observations. For our final results, 
we therefore adopt that NSs formed via AIC or MIC have short 
lifetimes (10 s years), and that NSs recycled in their original AIC 
systems are only mildly recycled, also with lifetimes of 10 years 
(case "2a" in [TJ. 

3.3 Recycling via stable MT 

A common understanding of recycled pulsar formation is that the 
NS is recycled through disk accret ion and has its magnet ic fields re- 
duced during the accretion phase dvan den Heuvell 1984) . Since the 
characteristic lifetime of a MSP is 10 11 — 10 13 years and is signif- 
icantly larger than a GC life-time (10 10 years), we do not consider 
in our simulations the evolution of MSPs towards a death zone after 
they are formed. 

The amount of material that needs to be accreted on to a NS 
to form a MSP is however not firmly established. In our simula- 
tions, we record for further analysis any event through which a NS 
acquires at least 0.01 Mq. We also distinguish mass gain through 
steady accretion from that which can result from a physical colli- 
sion with a red giant (DCE); a CE event during binary evolution; or 
a merger or physical collision with a MS star or a WD. 

The connection between stable mass-transfer and the bi- 
nary MSP production in the field seems to be reasonably well- 
understood and is reviewed bv lDelovd We outline only the 
main branches: 

• If the donor is a RG or an asymptotic-giant branch companion, 
then the resulting binary MSP has a wide orbital separation ("4b" 
at the Fig.[JJ>. 

• If the MT occurred in a NS-MS binary, then the evolution de- 
pends on the period at the start of the MT. Some systems will evolve 
towards shorter periods("3a" at the Fig. [TJ, and others to longer 
("4d" at the Fig. Q} (for a revi ew of the curr ent understanding of 
NS-MS LMXBs evolution see lDelovell2007l) . The progenitors to 
the short-period binary MSPs evolve slowly enough that they may 
be detected as (bright or quiescent) LMXBs ; systems with these 
orbital periods are ind eed seen in globulars i Heinke et alj|2005i3 ; 
lAltamirano et alj|2007l) . 

• If a binary system with a NS has evolved through a CE re- 
sulting in a significant shrinkage of a binary orbit, MT with a WD 
donor can start. This leads to the formation of an ultra-compact 
X-ray binary (UCXB) ("3b" at the Fig.[TJl. 

We do not discuss here the case of thermal-time scale MT that 
occurs with an intermediate-mass donor (2 ^ Md ^ 4Mq) in a 
close orbit at the start of the MT (0.3-2 days), as we find it un- 
likely for a NS in a GC environ ment to acquire an intermediate 
mass donor (c.f. lRasio et al.l200d) . We will return to this point later 
in §8 and only note here that the evolutionary life-time of an inter- 
mediate mass donor is not much longer than the cluster age when 
ECS NSs are formed. 

Even if the MT leads to spin-up, a NS does not appear as a 



radio MSP immediately even if has accumulated sufficient mass. 
While the MT rate remains high, systems continue to appear as 
LMXBs. In general, LMXBs with a likely MS donor do not show 
up as radio MSPs while th e companion mass is above ~ 0.2Mq 
(see Fig. 1 in |Pelovell2007h . 



3.4 Short-period radio bMSPs 

There are discrepancies between the observed orbital periods of 
short-period radio binary MSPs (bMSPs) with very low-mass com- 
panions and the modeled periods. Using a standard mass-radius 
relation for WDs during the MT evolution of a NS-WD binary 
gives periods several times smaller than those of observed bM- 
SPs with companions of similar masses. It has been shown that 
even the consideration of more realistic WD models does not 
change the result strongly - hotter donors don't, by themselves, 
prod uce the long periods needed to explain t he short-period bM- 
SPs dPelove & Bildstenll2003l ; iDelovd l2007h . If the mass-losing 
companion is a WD, some poorly understood process must be in- 
voked - e.g., strong tidal heating of the WD during MT (Rasio et al.l 
|2000|) . In that case, it is not possible to accurately predict the bMSP 
periods during the MT or at its end. 

There is no clear termination point for MT from a degenerate 
donor. As a result, the MT continues uninterrupted and within the 
Hubble time the companion mass can reach below 0.01 Mq. No 
observed radio bMSPs have been shown to have such low-mass 
companions, in globular clusters or the field. 

Several accretion-powered MSPs with P or b < 60 minutes 
have been id entified, consistent with (partly) degenerate He or C/O 
donors (e.g. Krimm et al. 2007), and several are known in globular 
clusters (e.g. iDieball et alJl2005D . However, no radio bMSPs have 
been observed with such short periods. These systems will spend 
Gyrs as transient LMXBs, with continually decreasing mass trans- 
fer rates ("3b" and "6a" in Fig.[T}. It is possible that the UCXBs 
in clusters remain in this configuration, never turning on as mil- 
lisecond pulsars, for a Hubble time. A sizable number of such sys- 
tems could remain hidden in globular clusters if their accretion out- 
burst intervals are sufficiently long. Alternatively, partly degener- 
ate (due to irradiation) donors will have a minimum companion 
mass for a given temperature, indicating that evaporation of the 
donor sta r and production of isolated millisecond p ulsars is a likely 
outcome dBildsterj2002l ; |Pelove & Bildsterj|2003l) . Other possibil- 
ities include very strong tidal heating of the companion leading to 
strong orbital expa nsion, producing the low-mass eclipsing pulsars 
dRasio et alJuOOOf) . or spin-down of the NS due to decreasing mass 
transfer rates (see discussion in |Pelovell2007l) . Or, current pulsar 
detection algorithms may not yet be sufficiently sensi tive to detect 
such short-period systems (e.g. iHessels et al.ll2007h . We do not 
know which of these fates befalls ultracompact LMXBs, although 
we know observationally that ultracompact radio bMSPs ("6a" in 
Fig. 1) have not yet been seen. We therefore count the numbers 
of ultracompact systems that have evolved past 55 minutes (the 
longest orbital per iod yet observed for an ultracompact LMXB, 
iKrimm et ai]|2007t) and tabulate them as ultracompact MSPs, but 
do not plot their final binary periods and masses; their true final 
fate is uncertain. 

In the case of MT in a NS-MS binary, if the resulting bi- 
nary has a short period, it will have a (partly) degenerate hydro- 
gen companion. Such a companion is larger than a hydrogen-poor 
WP, and better matches the observed periods of known low-mass 
eclipsing radio bMSPs than systems with WP donors. It is thought 
that the MT becomes unstable and leads to a merger of the two 
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stars when the hydrogen-rich companion bec omes less massiv e 
than ~ 0.016Mq (for analytic estimates see iKing etal.ll2005h . 
producing a single MSP ("4a" in fig. [D- However, current binary 
evolution codes lack the equation of state appropriate for low tem- 
peratures. Therefore, as of today, there are no direct mass transfer 
calculations that can proce ed below donor masses of 0.04 Msun 
(van der Sluvs et aTll2005al lbl), which is well above the masses of 
observed bMSPs with very low-mass companions. The main rea- 
son is that an equati on of state that is ap propriate for extremely 
low-mass stars (e.g., ISaumon et al.lll995l) is not yet incorporated 
in current binary evolutionary codes. 

Our code evolves a NS-MS binary down to very low masses 
assuming a simple mass-radius dependence for a hydrogen-rich de- 
generate donor. We find that this evolution ends as a merger when 
the donor mass is ~ O.O2A/0, where this value varies a bit depend- 
ing on the initial donor mass. While evolving, the system appears 
in a globular cluster as first a qLMXb and then as a bMSP for sub- 
stantial lengths of time. The values of the donor mass when the 
instability occurs, and the orbital periods during the qLMXB stage, 
could use refinement when improved evolutionary simulations are 
available. 



3.5 Recycling during fast ("dynamical") events 

We consider here mass gain which can result from a physical colli- 
sion with a red giant (DCE,case "4e", see Fig.QJ ; a CE event during 
binary evolution (case "4c"); or a merger or physical collision with 
a MS star or WD ("2b" at the Fig.Q]l. 

During a physical collision with a giant, the NS retains ^ 
O.OIMq of the giant's envelope dLombardi et alj|2006l) . This mat- 
ter has high specific angular momentum and the resulting accretion 
must proceed through an accretion disk, implying that a recycled 
pulsar can be formed. 

A CE event is also likely to lead to spin-up. As an example, it 
is believed that the double neutron star system PSR J0737-3039 i s 
a post-common envelope remnant l lDewi & van den Heuvef2004l) . 
where the first-formed NS has been spun up during a CE event to 
millisecond periods. 

It is less clear whether possible accretion as a result of a 
merger leads to recycling. Such mergers hav e been considered to 
produce only mild spin up dLvne et al]|l996l) . For a standard sce- 
nario we assume that physical collision of NSs with other stars 
leads only to the formation of short-living slow (isolated) pulsars. 



3.6 Post-recycling phase 

Once recycled, a single or binary MSP can participate again in dy- 
namical encounters, and form new binary systems whose proper- 
ties will depend on the type of binary formation: binary exchange 
("5a"), tidal capture ("5b") or another DCE ("5c"). This type of bi- 
nary millisecond pulsars, along with post-DCE "4e", should not be 
formed in the field population. Newly dynamically formed binaries 
with MSPs may evolve again through stable MT. 



4 GLOBULAR CLUSTER MODELS 

Our "standard" cluster model has an initial binary fraction of 100 
per cent. The distribution of initial binary periods is constant in 
the logarithm between contact and 10 7 d and the eccentricities are 
distributed thermally. We emphasize that about 2/3 of these bina- 
ries are soft initially (the initial binary fraction for hard binaries is 



about 20 per cent if the 1-D velocity dispersion is 10 km/s) and 
most very tight binaries are destroyed through evolutionary merg- 
ers. Our initial binary fraction is therefore comparable to the initial 
binary fractions that are usually adopted in iV-body simulations, 
10-15 per cent in hard binaries (soft binaries are short-lived in dense 
clusters but slow down the calculations considerably). For a more 
detailed discussion of the choice of primordial binary fraction, see 
llvanova et"al] d2005 j) . 

For single stars and prim aries we adopt the broken power law 
initial mass function (IMF) of lKroupal d2002l) and we assume a flat 
mass-ratio distribution for secondaries. The initial core mass is 5 
per cent of the total cluster mass and no initial mass segregation is 
assumed (an average object in the core is as massive as an average 
cluster star). At the age of 1 1 Gyr, which we adopt as a typical age 
of a Galactic GC, the mass of such a cluster in our simulations is 
~ 2.2 x 10 5 Mq for the metallicity of a typical metal-rich cluster 
(Z = 0.005) and is comparable to the mass of typical GCs in our 
Galaxy. In the case of a metal-poor cluster (Z=0.0005), the mass is 
~ 2.5 x 10° A/©. To make a comparison between cluster models 
of different metallicities we will show all results scaled to a cluster 
mass 2 x 10 5 Mq. As the metal-rich case observationally is more 
important - most bright LMXBs and MSPs detected so far are lo- 
cated there - we take for our standard model the value Z — 0.005. 

For our standard model we fix the core density at n c — 
10 5 pc~ 3 . The link to an observed luminosity density can be found 
using the number density-to-mass density ratio, which in our sim- 
ulations is ~ 2 Afg 1 at the ages of 7-14 Gyr, and mass-to-light 
ratio, w hich on average is ~ 2 .3 Mq/Lq and is 2 Mq/Lq for 
47 Tuc dPrvor & Mevlarill l993h. We adopt a half- mass relaxation 
time t r h = 1 Gyr l Harris! 19961) . The characteristic velocities are 
computed for a King model with dimensionless central potential 
Wo — 7 and the total model mass, giving a one-dimensional ve- 
locity dispersion o\ = 10 km/s and a central escape velocity 
«csc = 40 km/s. If, after an interaction or SN explosion, an ob- 
ject in the core acquires a velocity higher than the recoil velocity 
Wrec — 30 km/s, the object is moved from the core to the halo. The 
ejection velocity for objects in the halo is w e j,i, = 28 km/s. This 
standard model represents a typical dense globular cluster in our 
Galaxy. 

In addition to this "standard" model we also considered cluster 
models with the following modifications: 



a metal-poor cluster with Z 
different central densities 



0.0005 - "metal-poor"; 



pc 



n c = 10° pc~ J , n c = 10* pc" d 
"high-dens", "med-dens" and "low-dens" 



and n c = 10 3 
respectively; 

• a cluster with smaller velocity dispersion and smaller ejection 
velocity - "low-cr"; 

• shorter half-mass relaxation time "long-t r h"; 

• reduced initial binary fraction, — 50 per cent - "BF05" ; 

• magnetic braking prescription of Ra ppaport et al.l J 19831) - 
"fast MB"; 

• reduced efficiency of the common envelope, «ceA = 0.1 - 
"CE-reduced"; 

• natal kick distributio n is doubled peaked Maxwellian as in 
lArzoumanian et all d2002l) and NSs are produced only via core- 
collapse and double WD collisions (MIC), but not through EIC, 
AIC or evolutionary issued MIC - "oldkicks" 

• 47 Tuc-type cluster, charact erized by a higher den sity n c = 
10 5 4 pc" 3 (with M/L=2 as in dPrvor & Me\ 
sponds to p c — 10 4 ' 8 Lq pc -3 as inlHarrij 
ity Z = 0.0035, t lh = 3 Gyr dHarrislll996 



1996b - ) , metallic- 
11.5 km/s 



<7\ 
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ecc binary MSP 




UCXB& 
accretion 
powered 
binary MSP 
P<30 min 
KfGyr 



radio-quiet 
X-ray transient 
bMSP 
P> 30min 
few Gyr 



to n dtry / r& g y 



post-BE ecc binary MSP 
7-p> 10 Gyr, 
P = 1-10 3 days, ecc 
t -10 s -KfGyr 

coll 



post-TC binary MSP 
t p > 10 Gyr, 
P = 0.1-1 days 
t -5-10 Gyr 

mb J 







post-DCE ecc binary MSP 
7-p> 10 Gyr, 
P = 0.1-10 days, ecc 
t -ICF-Kfyr/r > Gyr 

gw coll 



Figure 1. The scenario for the formation and evolution of MSPs in globular clusters, rp is the life-time of a radio pulsar after formation, tmt i T mb i r gw , t co h 
indicate which time-scale determines the lifetime of the indicated binary: tmt - while MT lasts, T mb and r gw - until magnetic braking or gravitational 
radiation merges the binary, t co h - before the next dynamical encounter. The final fate of radio-quiet X-ray transient bMSPs ("6a") is not determined. 
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Table 1. Models. 





He 


'rh 


Z 




^CSC 


CO 
J b 


«ceA MB 


-^cluster 


Runs (number of stars) 


standard 


10 5 


1.0 


0.005 


10 


40 


100% 


1.0 IT03 


2 x 10 5 M© 


5 with 10 6 , 5 with 2 X 10 6 


metal-poor 






0.0005 












5 with 10 6 


high-dens 


10 6 
















5 with 10 6 


med-dens 


10 4 
















5 with 10 6 


low-dens 


10 3 
















5 with 10 6 


Iow-ct 








5 


20 








5 with 10 6 


long-t rh 




3.0 














5 with 10 6 


BF05 












50% 






5 with 10 6 


fast-MB 














RVJ 




5 with 10 6 


CE-reduced 














0.1 




5 with 10 6 


oldkicks 


















5 with 10 6 


47Tuc 


2.5 ■ 10 5 


3.0 


0.0035 


11.5 


57 






10 6 M Q 


5 with 2 X 10 6 


Terzan 5 


8 ■ 10 5 


0.93 


0.013 


11.6 


49 






3.7 x 1O 5 M 


5 with 1.7 X 10 6 



Notations: n c - core number density (pc~ 3 ); t T ^ - half-mass relaxation time (Gyr); Z - metallicity; a\ - one-dimensional velocity dispersion (km/s); v csc - 
central escape velocity (km/s); fP - initial binar y fraction (see text for co rresponding bina ry fraction of hard bina ries); oqeA - CE efficiency; MB - adopted 

prescription for magnetic braking (IT03 is from llvanova & Taaml 120031) and RVJ is from lRappaport et alj ll983h); M c i uster - all th e results for this model 
are scaled to this adopted cluster mass. "Oldkicks" model uses the distribution of natal kick velocities from lArzoumanian et alj2002l) . Symbol "■" means that 

the value is the same as in the standard model. 



JPrvor & Mevlarl 19931) . v BSC = 57 lWebbinkll985h . with the recoil 
velocity of 51 km/s and u e j,h = 32 km/s (King model Wo = 12)) 
- "47 Tuc"; 



Table 2. Production of NSs, field population. 



Terzan 5 cluster, n c = 10 pc 
s in Heinke et all 120031)). metallic ity Z 
0.9 Gvr lHarris|[T996h ax 



(fie 



10 J 



0.013 lOriglia & RicH 
120041). Uh = a 9~Gvr lHarrisll 19961) o\ = 11.6 km/s, v esc = 49.4 
lWebbinkll985r) . with the recoil velocity of 39 km/s and u e j,h = 29 
km/s (King model W = 8.3) - "Terzan 5"0 

For the complete list of models see Table[T] 

We expect that only a few LMXBs can be formed in every 
cluster with a typical mass of 200,000 Mq through 1 1 Gyr of dy- 
namical evolution, and therefore we perform 5 runs for each model 
with 10 6 initial stars and for the standard models - five runs with 
2 x 10 6 stars, to check that statistics for rare objects like LMXBs 
are not different when the resolution is increased (for most of our 
results presented below we did not find statistical differences be- 
tween high and low resolution models, and fluctuations for the rate 
events are smaller in the case of high resolution). For Terzan 5 and 
47 Tuc models, though, we do all runs at high resolution. In partic- 
ular, in the case of Terzan 5 1.7 x 10 6 stars provide at 1 1 Gyr the 
mass comparable to the estimated observed mass of 370,000 Mq. 
In the case of 47 Tuc the reason for a higher resolution is that the 
relative core mass is smaller than in the case of a standard cluster, 
due to longer t r h. 



5 NEUTRON STAR FORMATION AND RETENTION 
5.1 Formation 

To estimate how many NSs can be formed and retained in GCs, 
we first studied stellar populations without dynamics, considering 
separately populations with only primordial single stars and only 



3 The recent study of lOrtolani et alj 120071) suggests a shorter distance, and 
thus higher central density for Terzan 5 (p c = 10 5 ' 61 L©/pc 3 ). If this study 
is correct, then Terzan 5's characteristics may be more closely approximated 
by our "high-density" cluster. 
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Z=0.001 
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Z=0.02 
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400 






1.2 
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424696 


oldkicks 
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19.3 


21.9 


339669 


binary 
















Z=0.0005 


3079 


545 


59 


14 


19 


329 


249775 


Z=0.001 


3056 


553 


60 


15 


15 


332 


242601 


Z=0.005 


2750 


576 


58 


16 


14 


335 


219236 


Z=0.02 


2463 


406 


33 


20 


12 


240 


180582 


oldkicks 


2929 








58 


58 
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Notations for channels - CC - core-collapse supernova; EIC - evolution 

induced collapse via electron-capture supernova; AIC - accretion induced 
collapse; MIC - merger induced collapse. Numbers are scaled per 200,000 
Mq stellar population mass at the age of 11 Gyr. 7V|q and are 
the retained numbers of NSs formed via core-collapse and through all 
channels, if the escape velocity is 40 km/s. /Vns are the total numbers of 
formed NSs in simulated populations. 



primordial binaries. We considered several metallicities - from so- 
lar (as in the Galactic field) to that of a typical metal-poor clus- 
ter: Z=0.02, 0.005, 0.001, 0.0005. In addition, for comparison with 
studies done by different researchers in the past, we examined a 
model (with Z=0.005) with the previously widely accepted natal 
kick distributio n, 2 Maxwellians with a lo wer peak for kick veloc- 
ities at 90 km/s lArzoumanian et aT]|2002h . In this model we adopt 
that NSs are produced only via core-collapse or through collisional 
WD mergers, no evolutionary ECS. The results for NS production 
via different channels are shown in Table|2] 

The production of NSs via core-collapse SNe (CC NSs), per 
unit of total stellar mass at 11 Gyr, decreases as metallicity in- 
creases. The difference between the Galactic field case and a metal- 
poor case is about 20 per cent, while the difference between metal- 
rich and metal-poor clusters (Z=0.005 and Z=0.0005) is about 10- 
15 per cent, for primordial single or binary populations. 



10 N Ivanova et al. 




v [km/s] 

escL / J 



Figure 2. The retention fractions as a function of escape velocity (for stel- 
lar evolution unaffected by dynamics). Dotted and dash-dotted lines show 
the retention fractions for single and binary populations, core-collapse NSs 
only. Solid and short-dashed lines show the total retention fractions for sin- 
gle and binary populations, all NSs. Long-dashed line shows the total reten- 
tion fraction of a binary population with the old kick distribution. 

EIC in single stars for metal-rich populations occurs in stars 
of higher masses. The range of initial masses is, however, the same 
for both populations and is about O.GMq. As a result, the number 
of NSs produced via EIC (EIC NSs) from the population of single 
stars is smaller in the case of a Galactic field population than in 
the case of a GC population by 30 per cent, while the difference 
between metal-rich and metal-poor clusters is only 5 per cent. This 
is in complete agreement with the adopted power law for the IMF. 
In binaries, the mass interval where EIC could occur is expanded 
compared to single stars of the same metallicity. 

NSs produced via AIC and MIC also have similar production 
rates in metal-rich and metal-poor populations. We note that MIC 
gives the smallest contribution to the production of NSs in popula- 
tions evolved without dynamics. 

It can be seen that once we adopt that a NS can be created 
via an ECS, the number of formed NSs via CC is depleted - some 
stars would be evolving via core-collapse, but have experienced the 
condition for a EIC before that. 



5.2 Retention 

In Fig. [2] we show the retention fraction of NSs at formation as a 
function of fixed escape velocity the whole population, consider- 
ing populations of metal-poor primordial singles and binaries; CC 
NSs and the whole NS population are shown separately. It can be 
seen that, with the adopted natal kick distribution, CC NSs from 
a single star population are rarely retained (see also Tabled- In a 
population of primordial binaries, the fraction of retained CC NSs 
is higher than for single stars. None the less, the resulting number 
of retained CC NSs is significantly smaller than we get from the 
other NS formation channels. We conclude that if indeed EIC, AIC 
or MIC are not accompanied by a large natal kick, the total number 



of retained NSs is high and only a negligible number of CC NSs 
is present. As the EIC/ AIC NS masses are only 1.26A/©, this may 
lead to a NS mass function skewed towards low values. 

With the old distribution for natal kick velocities, reasonable 
numbers of NSs can be produced and retained in a typical GC, 
while a massive cluster like 47 Tuc could retain as many as 600 
NSs (see al so Fig.[2]>. This res ult agrees with the retention fractions 
obtained in lPfahl let al] j2002h . 



5.3 Effects of dynamics 

The effects of dynamics on the production and retention of NSs are 
illustrated in Table|3] The number of formed CC and EIC NSs is in 
agreement with the retained fraction of the stellar population where 
binaries have been partially depleted and some of the SN progeni- 
tors were located in the halo, where the escape speed is lower than 
in the core. We find that the production of NSs via AIC and MIC 
is strongly enhanced by dynamical encounters compared to a field 
population. As a result, those formation channels together become 
comparable to the EIC channel, and produce many more retained 
NSs than the CC channel. A metal-poor cluster, following the trend 
of field populations, possesses a few per cent more NSs. 

As noted above (§3.1), MIC could lead to the formation of 
magnetars, which also appear as anomalous X-ray pulsars (AXPs). 
Currently, there are no AXPs detected in old GCs, and only 
one A XP is detected in a young "super star cluster," Wester- 
lund 1 dMuno et alfcOOfj) . On the other hand, all known AXPs are 
younger than 10 5 yr. In this case, AXPs could be found in GCs 
only if their formation rate is above 1 per 200 (GCs) per 10° yr, or 
at least 500 per typical cluster, over the cluster lifetime. Our forma- 
tion rates are significantly smaller and therefore the hypothesis of a 
MIC connection to AXPs does not contradict our results. 



5.4 Location 

We consider the position of retained NSs in the cluster. The progen- 
itors of CC NSs are mainly massive stars with masses > 20M©. 
However a significant fraction of progenitors of EIC NSs, espe- 
cially those that evolved from a primordial binary, have intermedi- 
ate initial masses — as low as 4 Mq. The evolutionary time-scales 
for NSs to be formed are ~ 10 7 yr for an CC NS and ~ 10 s yr 
for an EIC NS. These are only ~ 2 per cent and ~ 10 per cent, 
correspondingly, of the cluster half-mass relaxation time i r h, where 
t r h = 10 9 yr in a typical cluster. Even massive stars do not experi- 
ence strong mass segregation during ~ 0.02i r h, the same is true for 
sta rs of intermediate m asses that produce EIC NSs (see, e.g., Fig.7 
in iGurkan et al.l|2004[ , , where the mass segregation with time in a 
cluster with stars between 0.2 Mq and 120 Mq is shown). When a 
NS is formed, the mass of the star is significantly decreased. As a 
result, a significant fraction of NSs is formed in the halo, and many 
will remain there. As many as ~ 60 per cent of all EIC NSs in the 
case of a typical cluster (t r h = 10 9 yr) remain in the halo at 11 
Gyr. AIC and MIC NSs are also present in the halo, but in a lesser 
proportion, as they are mainly produced in a dense environment. 
Many of the halo AIC and MIC NSs have recoiled from the core as 
a result of dynamical encounters. A significant difference with the 
distribution of NSs over the clusters can only be expected for ini- 
tial mass segregation, or if the initial trh was much shorter than the 
current resultin g in faster mass segre gation at early stages. Ob- 
servations find (e.g. lGrindlav et al.l2002h that the radial distribution 
of MSPs within globular clusters in 47 Tuc is consistent with the 
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Table 3. Retained NSs and their locations at 1 1 Gyr, for different globular cluster models. 



model 




Core 






Halo 






CC 


EIC 


AIC 


MIC 


CC 


EIC 


AIC 


MIC 


standard 


3.5 ± 1.5 


65.6 ±9.3 


21.1 ±3.7 


10.3 ±4.0 


4.8 ± 1.7 


91.4 ±11 


18.0 ±3.9 


9.7± 1.7 


metal-poor 


3.6 ± 1.7 


71.5 ±12 


30.4 ±3.2 


8.4 ±4.5 


2.3 ± 1.1 


91.7 ±9.2 


20.9 ±4.8 


10.9 ± 3.8 


high-den 


3.3 ± 1.7 


65.5 ±7.8 


42.7 ±5.8 


20.7 ±6.1 


4.2 ±3.3 


111 ± 12 


29.0 ±5.2 


12.8 ± 3.4 


med-den 


3.3 ±2.1 


63.7 ±5.2 


18.6 ±4.5 


6.4 ± 3.0 


4.8 ± 1.5 


85.6 ±8.0 


14.3 ±3.1 


9.4 ± 3.7 


low-den 


3.2 ± 1.3 


65.1 ± 5.9 


14.2 ±3.3 


9.8 ±0.9 


3.8 ± 1.3 


91.2 ± 11 


15.4 ±2.1 


7.0 ±2.0 


lowsig 


0-5±8;| 


15.5 ±3.2 


5.5 ±2.3 


10.7 ±2.5 


9+ 11 


19.4 ±4.2 


4.3 ±1.7 


11.5 ± 1.8 


longtrh 


1.4 ±1.0 


30.0 ±4.7 


9.3 ±2.1 


3.0 ± 1.2 


5.3 ±0.9 


124 ± 15 


22.5 ±5.5 


15.5 ± 3.3 


BF05 


1.8 ± 1.2 


64.1 ± 5.0 


15.0 ±3.1 


6.7 ± 2.5 


3.9 ± 1.0 


88.6 ±6.7 


11.4 ±2.7 


6.2 ± 1.7 


fast-MB 


3.8 ± 1.6 


77.3 ±6.0 


22.9 ±3.3 


10.6 ±4.3 


5.1 ± 1.9 


93.7 ±6.9 


17.2 ±3.3 


8.3 ±3.2 


CE-reduced 


2.9 ± 1.9 


72.2 ±5.6 


12.7 ±2.8 


9.4 ± 1.9 


3.6 ±2.6 


99.2 ±4.0 


4.5 ±1.7 


11.7 ± 1.9 


oldkicks 


14.6 ±3.9 


0.0 ±0.0 


0.0 ±0.0 


u.oze 3 


21.3 ±4.0 


0.0 ±0.0 


0.0 ±0.0 


0.0 ±0.0 


Tuc-47 


11.8 ± 7.5 


215 ± 20 


91.8 ± 23 


27.7 ±8.0 


38.3 ± 7.9 


598 ± 41 


102 ±4.9 


69.3 ± 12 


Terzan-5 


5.1 ± 1.4 


132 ±8.9 


91.9 ± 11 


47.8 ±4.6 


8.9 ± 1.4 


143 ± 10 


38.0 ±6.6 


22.7 ±4.5 



Notations for channels as in Table|2] For all GC models, except 47 Tuc and Terzan 5, the numbers are scaled per 200 000 Mq stellar population mass at the 
age of 1 1 Gyr; for 47 Tuc the numbers are given per its total mass taken as 10 6 Mq , for Terzan 5 - per 370 000 Mq . Numbers show the results averaged 

over all runs for each cluster model (for the number of runs see Table[T}. 



produ ction of all MSPs within the core (but cf. iGurkan & Rasid 
l2003h . 

6 FORMATION OF BINARIES WITH NEUTRON STARS 
6.1 Primordial binaries with NSs 

The probability for a NS to remain in a primordial binary depends 
strongly on its formation channel (see Table[4]l. By definition, MIC 
leads to the formation of a single NS. AIC occurs in a very tight 
binary. Accordingly, a NS formed this way will most likely remain 
in a binary. Among CC NSs or EIC NSs, only several per cent re- 
main in binaries immediately after NS formation. Binaries that sur- 
vive SN generally have relatively massive secondaries that are also 
likely to experience SN, or the binary may evolve through dynam- 
ically unstable mass transfer resulting in a merger, with stronger 
binary depletion in the case of initially more massive binaries with 
a CC NS. As a result, very few NS binaries formed through those 
channels remain in binaries. Post-AIC binaries are the GC primor- 
dial binaries most likely to contain a NS. The low probability for 
CC NSs and EIC NSs to be members of a primordial binary does 
not change strongly with metallicity. 

The formation rate of AIC NSs in GCs shows that as many 
as 40 NS-MS binaries can be formed in a typical cluster. Most of 
them will evolve through the LMXB stage before the age of 1 1 Gyr, 
producing as a result 20-30 recycled pulsars, although those MSPs 
may have short lifetimes and not be observable today (see ij3,2| >. 
Overall, NSs that remain in binaries may evolve through the mass- 
transfer phase, but most such events occur very early in the cluster 
history, mainly within the first 2 Gyr. In more detail, per sample of 
270 000 NSs formed in primordial binaries (Z=0.005), no binaries 
started the MT with a MS star after the stellar population reached 
the age of 8 Gyr; for ages between 10 and 1 1 Gyr we encountered 
only 1 LMXB with a RG and 2 LMXBs with a WD companion. 
From these stellar population runs it is clear that it is highly un- 
likely that LMXBs observed in GCs of our Galaxy today are pri- 
mordial. Exceptions to this can occur only if a primordial binary 
had participated in an encounter, retaining both its companions but 
changing binary separation and/or eccentricity, resulting in mass 
transfer at older cluster ages. However, simply being in a binary 



Table 4. Binary fractions for NSs in primordial binaries. 



z 


after NS formation 




at 1 1 Gyr 






CC 


EIC 


AIC 


CC 


EIC 


AIC 


0.0005 


2.9% 


3.0% 


92.3% 


0.28% 


1.38% 


78.6% 


0.001 


2.8% 


3.4% 


93.6% 


0.25% 


1.41% 


81.8% 


0.005 


2.6% 


3.4% 


96.0% 


0.18% 


1.54% 


88.0% 


0.02 


1.4% 


4.2% 


95.2% 


0.12% 


1.49% 


85.9% 



The fractions of NSs that remain in primordial binaries immediately after 
NS formation (columns 2 v 4) and at 1 1 Gyr (columns 5 -f- 7), for each 
formation channel. 

immediately after NS formation strongly increases the chance for a 
NS to participate in an encounter and, accordingly, to get another, 
non-primordial, companion. The fraction of NS binaries with a NS 
formed via AIC is therefore expected to be higher than for other 
NS channels. 

The number of double neutron star (DNS) systems formed in 
the field is very low, typically just a few DNSs per 200,000 Mq at 
1 1 Gyr. For a typical GC this number will be smaller since some 
of the primordial binaries will be destroyed dynamically and most 
DNSs will be ejected. In fact, more than 99 per cent of the first- 
born NSs in primordial DNS systems are CC NSs, and, accord- 
ingly, have high natal kicks. The formation of primordial DNSs 
occurs very early in the GC evolution. Most of them are very tight 
and merge quickly, leaving very little chance for one to be observed 
today. As this paper only uses the results of the NS formation and 
evolution in the field as a reference point, we refer the reader to the 
complete study of DNS formation channels from primordial bina- 
ries, rates, appearance with time, metallicity dependence and other 
details provided in Belczynski et al. 2007 (in prep). We conclude 
that essentially no primordial DNSs can be present in Galactic GCs 
at the current time. 



6.2 Physical collisions 

A binary with a NS and a WD can be formed vi a physical coll i- 
sion with a g i ant (s ee £|2.2.U for earlier works see lVerbuntl dl987l) ; 
iDavies et al This binary formation mainly works during 
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Table 5. NS binaries formed in scattering experiments. 





Z 


17 


DCE 




TC 






km/s 


Tot 


MT 


Tot 


MT 


standard 


0.005 


10 


3.68% 


1.75% 


2.35% 


1.20% 


metal-poor 


0.0005 


10 


3.25% 


1.30% 


1.93% 


0.95% 


low-cr 


0.05 


5 


6.80% 


3.61% 


9.24% 


4.37% 



The table shows the fraction of NSs that successfully formed a binary via 
physical collision with red giants (DCE) or tidal capture (TC). Tot is the 
total number and MT is the fraction of NSs that not only formed a binary, 
but also started mass transfer before 1 1 Gyr. 



100 



10 



0.1 



0.01 




0.2 



0.4 



0.6 



Figure 3. Companions in binary systems formed via physical collision with 
a red giant during the scattering experiment (Z=0.001). Triangles are WDs, 
circles are He stars; solid symbols denote the systems that started the mass 
tranfer before 1 1 Gyr. 



single-single star collisions, although it may occur during a binary- 
single encounter as well. As we showed above, the number of NSs 
that can be retained in globular clusters is about a couple hundred 
for a typical globular cluster. The cross-section for two stars that 
have a relative velocity at infinity Voo to pass within a distance r m ax 
is 



1 + 



2G(mi +m 2 ] 



), 



(9) 



' max t^oo 



where the second term is gravitational focussing and in, stellar sys- 
tems with low velocity dispersion like GCs, it exceeds the first term. 
The time-scale for a NS to undergo an encounter with a red giant 
can be estimated as r co u = V /rg nav^ , where /rg is the relative 
fraction of red giants in the stellar population: 



Tcoll 



1.6 x 10 11 yr vion B (Mns + A/rg) 1 Rml x fnG > 



(10) 



where A/ns and Mrg are the masses of the NS and red giant in 
Mq, and -R m ax = »" m ax/R© is the maximum distance between 
the stars during the en counter that leads to the binary formation. 
lLombardi et al .1 l |2006i) shown that r max ~ 1.5-Rrg, where Rug 
is the radius of the red giant. Averaged over time, Rn.G is a few 



solar radii dlvanova et ai]|2005dh , and /rg ~ 5 per cent. Using 
eq. d 10b we estimate that during 1 1 Gyr about 5 per cent of NSs can 
participate in the formation of binaries through physical collisions. 

As the estimate shows, having about 200 NSs in our mod- 
els of GCs (and about 400 for the models with the highest reso- 
lution, 2 x 10 6 stars), we can form a maximum of 10 binary sys- 
tems through the whole cluster simulation. To study statistically 
the characteristics of post-encounter binaries, we performed a scat- 
tering experiment; we considered a population of 10 6 single stars, 
with an IMF from O.2M to 3.0M© and 10000 NSs of 1.47Y . All 
stars are initially placed in the core with n c = 10 5 pc -3 . We con- 
sidered three different models - standard, low metallicity, and with 
low velocity dispersion (see also Table[5]l. In these simulations only 
single-single encounters were allowed. 

The scattering experiment shows that 3.7 per cent of all NSs 
successfully formed a binary via physical collision with a RG or 
an AGB star. This fraction slightly decreases when the metal-poor 
population is considered and increases by a factor of about 2 (in 
accordance with eq. 1 10b in the case of lower velocity dispersion. 
The fraction of NS-WD systems that started MT before 1 1 Gyr is 
high, ~ 50 per cent (see Tabled- 

In Fig.[3]we show the masses of companions in post-collision 
binaries vs. orbital periods. We form about the same number of sys- 
tems with WD companions and stripped He star companions (from 
collisions with AGB stars). Systems with He star companions are 
mainly formed within the first 2 Gyr and start MT mostly within 3 
Gyr. After 3 Gyr both the formation rate and the appearance rate 
are slightly higher for systems with WD companions. The evolu- 
tion of mass-transferring systems with He star companions is a bit 
different from those with WD companions. They are a few times 
wider (for the same companion mass) and do not live as long in the 
mass-transferring phase at high luminosities. 

In our numerical simulation of different models of GCs the 
actual rates are smaller. For instance, not all NSs are immediately 
present in the core at the time of their formation, and about half of 
them are still in the halo at the age of 1 1 Gyr. As a result, our "stan- 
dard" model shows about 2 per cent formation rate per NS present 
in the core (or 1 per cent per all NSs in the cluster), producing only 
2 systems. The scatter between several realizations is big. In some 
simulations we formed only half as many binaries via this channel, 
and, in one simulation, 3 times more binary systems were formed. 
For our standard model, no simulations failed to produce binaries. 
We note that the low-cr cluster model, unlike the scattering results, 
does not produce twice as many NS-WD binaries as the standard 
model - as the low-cr cluster retains less NSs than the standard. 



6.3 Tidal captures 

In our previous study, devoted to t he formation of catac lysmic vari- 
ables (CVs) in globular clusters dlvanova etai]|2006h . we found 
that tidal capture plays a relatively small role in the formation of 
dynamical binaries, only ~ 1 — 2 per cent of the total formation 
rate. Two main reasons were responsible: (a) a significant fraction 
of CVs have been formed directly from primordial binaries and (b) 
both physical collisions and exchange encounters occurring with 
WDs form tighter binaries, increasing their chance to start the MT. 
For the case of tidal captures with NSs, the results are different. 

In our scattering experiments, we find that only ~ 2 per cent 
of NSs form a binary via TC (see Table |5); with a small decrease 
in the case of a metal-poor population and an increase (4 times) 
for the lower velocity dispersion case. The latter effect is due to 



»"max being a function of the re lative velocity, and in creasing as the 
relative velocity decreases (see llvanova et ail d2006l) ). 

The rate of tidal captures is flat both with time and companion 
mass. In the case of the metal-poor population, most TC systems 
that start MT before 1 1 Gyr have formed before 5 Gyr. For our stan- 
dard metal-rich population, among systems formed in the past few 
Gyrs, only those with MS companions > 0.85A4© are able to start 
MT. We conclude that < 1 per cent of NSs can form a LMXB via 
tidal captures, and that the number of LMXBs mainly depends not 
on the current cluster properties, but on the dynamical conditions in 
the cluster core more than 5 Gyr ago. For metal-rich clusters, there 
is an exception for initially more massive MS companions. 

In our numerical simulations the formation rate is even 
smaller, due to the same reasons as in tj6.2l An average probability 
in our "standard" model for a NS to participate in a TC is 0.6 per 
cent, and half that to form an LMXB via this channel. None the 
less, even being very small, the TC formation channel for LMXBs 
is comparable to that of physical collisions (see for more details 
0. 
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Table 7. Average number of LMXBs with different donors 



model MS RG WD 



standard 


1 


85 ± 0.56 





19+0.16 





46 ± 0.41 


metal-poor 


1 


22 ±1.05 





10 ±0.07 





16 ± 0.16 


high-den 


7 


31 ±2.21 





on+0.27 
zu=c 0.20 


6 


38 ± 3.60 


med-den 
low-den 
low-cr 





3 


46 ± 0.40 

ci +0.63 
Oi± 0.51 

13 ± 1.25 







no+0.04 
01 ±0.01 

07 ±0.06 







no^o.12 

U8± 0.08 
J " L:E 0.11 

12 ± 0.11 


long-t rh 
BF05 



1 


84 ± 0.73 
12 ±0.32 






Oq+0.18 
ua:c 0.09 

08± o:o8 






11+0.16 
11=c 0.11 

44 ± 0.33 


fast-MB 


1 


38 ±0.93 





1n +0.11 


1 


09 ± 0.75 


CE-reduced 


1 


98 ±1.33 





06 ± 0.06 





23 ± 0.16 


oldkicks 





nq+0.49 
z9:t 0.29 














47Tuc 


9 


28 ±4.10 





13 ±0.09 


i 


57 ±3.03 


Terzan 5 


7 


16 ±1.95 





14 ±0.11 


7 


73 ± 3.03 



Average number of LMXBs present at ages 11 ± 1.5 Gyr, for different 
donors (RG includes also subgiants and Hertzsprung gap stars). 



6.4 Exchange encounters 

As can be seen from eq.[l0] due to a larger r ma x, binary encounters 
are much more likely than single star collisions. However, com- 
pared to the case of physical collisions and TC, where very tight 
binaries are formed, a smaller fraction of the successful exchange 
encounters are expected to result in a binary with MT. The rea- 
son is that when the binary exchange occurs, the binary separation 
in the formed post-exchange binary roughly scales with the pre- 
exchange binary separation as the ratio of the new compani on mass 
to the mass of the replaced companion l lHeggie et al]|l99ot) . There- 
fore, when a single NS becomes a member of a binary system at 
the age of several Gyrs, the post-exchange system expands, as the 
mass of the replaced companion is always less than the mass of 
the NS. In addition, when a binary encounter involves a very tight 
binary, which would be likely to start MT soon even in the case of 
post-exchange expansion, a physical collision during the binary en - 
counter is more likely than a simple exchange dFregeau et al.l20Q4h . 

Overall, in our standard model ~ 20 — 25 per cent of NSs 
in a cluster (40-50 per cent if we consider only NSs in the core) 
will successfully form binaries through exchange encounters. ~ 80 
per cent of post-exchange binaries will have a MS companion, and 
~ 15 per cent a WD companion. Only 1.5 per cent of NSs will 
begin MT with a MS companion, and only 0.2 per cent with a 
WD companion. The role of binary exchanges in the formation of 
NS-WD LMXBs is therefore negligible compared to physical col- 
lisions. 



7 LOW-MASS X-RAY BINARIES 

From the discussion above, one can see that the expected forma- 
tion rate of LMXBs is rather small. For all dynamical formation 
channels that involve a NS directly, in a typical dense metal-rich 
cluster, the combined rate is only ~ 3 per cent per NS. This is 6 
LMXBs per 11 Gyr of the cluster life, where 2 LMXBs are with 
a WD companion and 4 LMXBs are with a MS or a RG compan- 
ion. A small number of LMXBs may come from primordial, but 
dynamically modified binaries (where the encounters changed the 
binary eccentricity or separation). In addition, however, some frac- 
tion of LMXBs is expected t o be provided by dyn amically formed 
binaries with heavy WDs (see llvanova e t al. (2006)). Such binaries, 



if evolved through AIC, create binaries that reinstate MT on a NS 
soon after AIC, with a MS, WD or RG donor. From our current 
simulations, we find that the channel of LMXB production from 
AIC of dynamically formed WD binaries produces more LMXBs 
than any of the other dynamical LMXB production channels dis- 
cussed above. Overall, in a typical cluster it is 2-3 times more ef- 
ficient during 9.5-12.5 Gyr than physical collisions, tidal captures 
and binary exchanges (see Table [6}. The exceptions are only very 
high density clusters. We also note that if the majority of LMXBs 
are produced through AIC, but these binaries do not produce long- 
lived MSPs ( fl3~2t , then the production rates of LMXBs and MSPs 
may not be closely coupled. However, LMXBs produced by AIC 
may also have shorter lifetimes due to the exhaustion of the com- 
panion's mass during the pre- AIC MT stage. 

In Fig. [4] we show when LMXBs (with different compan- 
ions) appear during the cluster evolution. We show only the cluster 
ages after 4 Gyr, as before that time mainly primordial LMXBs 
appear. The connection between the production rates of LMXBs 
and the numbers that are present at a particular time in a GC is 
provided by the life-time tlmxb of LMXBs with various donors. 
An average tlmxb for NS-MS LMXBs is about 1 Gyr. Depend- 
ing on the metallicity and the initial donor mass, a system can 
be persistent 5-40 percent of the MT time for metal-rich donors 
> 0.6M© and transient during the rest; for donors of lower metal- 
liciti es or smaller m asses, a NS-MS LMXB will be transient all the 
time jlvanovdl2006T) . and therefore more likely seen as a qLMXB 
rather than as a bright LMXB. An LMXB with a red giant com- 
panion or a companion in the Hertzsprung gap is very short-lived, 
10 5 — 10 7 yr, and in only very rare cases can they live as long as 
10 8 yr. In the case of NS-WD LMXBs - ultra-compact X-ray bi- 
naries (UCXBs) - the total tlmxb is a few Gyr, however the time 
during which the system is persistent and has an X-ray luminosity 
above 10 36 erg/s is only ~ 10 8 yr. UCXBs at ~ 10 8 years will 
fall below Lt, i = 1.3 x 10 36 ergs /s. According to the theory of 
irradiated disks jPubus etal.ll 19991) and the evolution of UCXBs 
dPelove & Bildstenll2003Tr 10 36 er gs/s is roughly the crit ical Lx 
for maintaining persistent accretion, lin't Zand et al .1 J2007I) present 
a table of 40 persistent X-ray bursters, which shows a striking low- 
L x pileup and cutoff at Lx = 1.3 x 10 36 ergs/s, or 0.5% of the 
hydrogen-poor Lsdd (incorporating their relativistic correction). 
Incompleteness in Galactic surveys does not fully explain this cut- 
off, as the PCA bulge scans monitor many sources at half this flux. 
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Table 6. Average appearance rate of LMXBs formed via different dynamical channels. 



model 


AIC 


TC 


DCE 


BE/MS 


BEAVD 


standard 


1.49 ±0.54 


0.07±°;i? 
0.11±°;}f 


0.14±g;a 


0.51±°;ff 


0-06±g;l° 


metal-poor 


2.40 ±0.54 




0.67±g; 6 73 7 


0.0 


high-den 


3.62 ±1.42 


0.97 ±0.76 


1.27 ± 0.66 


1.81 ±0.71 


0.0 


med-den 


0.54 ±0.19 


11+ 015 


0.0 


o.ii±g:» 
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Number of LMXBs that appear per Gyr at age 11 ± 1.5 Gyr, for different formation channels: AIC in a dynamically formed WD-binary, TC - tidal capture, 
DCE - a dynamical common envelope event (physical collision with a red giant), BE/MS - binary exchanges when a MS companion is acquired, BE/WD - 

binary exchanges when a WD companion is acquired. 
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Figure 4. Number of appearing LMXBs per Gyr with different donor types (with all donors, with MS donors, with WD donors and with donors that are RGs, 
subgiants or are in the Hertzsprang Gap). The left panel shows the results for the standard model (solid line) and the metal-poor model (dotted line). The right 
panel shows the results for 47 Tuc (solid line) and Terzan 5 (dotted line). 



Therefore we adopt Lx = 1-3 x 10 36 ergs/s as the break between 
persistent and transient behavior for UCXBs. 

Combining the number of appearing LMXBs with their life- 
times in the simulations, we find that a typical dense cluster at 
11 ± 1.5 Gyr age can contain up to 2 LMXBs with a MS com- 
panion and up to 1 LMXB with a WD companion (see also Ta- 
ble [7] note the very large scatter). We specify, that Table [7j shows 
the average expected number of LMXBs - both bright LMXBs and 
qLMXBs. We choose a cutoff period of 55 minutes for identifying 
UCXBs as LMXBs. We justify this by noting that th e longest pe- 
riod U CXB yet identified has a period of 55 minutes dKrimm et al.l 
|2007|) . and that the evolution of UCXBs will produce very large 
numbers of longer-period systems which have not yet been seen. 
As expected, the number of LMXBs shows a strong dependence 



on the core density; dependence on the velocity dispersion and the 
half-mass relaxation time is also observed. 

A strong dependence of LMXB production on cluster density 



Heinke et alj 


2003) 


Jordan et alj 


2004 



(e.g. JVerbunt & Hut|[l987l ; IPoolev et alj|2003l : 

and in extragalactic globular clusters (e.g., 
ISivakoff et alj [2007h . The interpretation has been that TVlmxb 
scales roughly with the encounter frequen cy T — p^.r^/a, where 
r c is the core radius dVerbunt & Hutl 19871) . 

The observations of LMXBs in our galaxy suffer from low 
number statistics-only 12 clusters have LMXBs-so attempts to un- 
derstand the LMXB formation rate from them encounter difficul- 
ties (e.g. iBregman et al]|2006l) . Observations of qLMXBs in our 
galaxy also suffer from low number statistics, as only ~two dozen 



Binaries with NS in globular clusters 15 




Figure 5. The specific collision numbers 7 and numbers of LMXBs in sim- 
ulated clusters. The solid line corresponds to tilmxb = (l±l) + O.O287. 
The error bars correspond to the scatter in our simulations. 



clusters have been observed with Chan dra sufficiently to iden- 
tify most qLMXBs. Ipoolev & HuJ J2006h and lHeinke etall d2006t) 
study overlapping samples of clusters using different assumptions. 
Their estimates of the density dependence of qLMXB production 
for the case of no metallicity dependence also overlap, suggesting 
a de nsity dependence lar ger than suggested by T, but with large er- 
rors. |Hejnke_etaJ 12006J) also investigate (inconclusively) the pos- 
sibility of metallicity dependence, which decreases the required 
density dependence. 

Observations of LMXBs in globular clusters of other galax- 
ies do not suffer from low statistics, and thus have clearly demon- 
strated a strong dependen ce of LMXB production on metal licity 



( Kun du et al 



Sivakoff ct al 



20021, |2007h . Analyses by Ijordan et alj d2004l) and 
d2007l) both find that LMXB production has a den- 
sity dependence lower than T. However, it is extremely difficult to 
determine the structural properties of extragalactic globular clus- 
ters with high accuracy, so uncertainty in the density dependence 
remains. 

In Fig.[5]we demonstrate how the specific number of LMXBs 
hlmxb = -/Vlmxb /M depends on the speci fic collision frequenc y 
7 = r/M, taken per units of 10 6 M , as in lPoolev & Hutl d2006l) . 
To calculate T in the same way as observers, we extracted values 
of p c and r c from our models at the cluster age of 11 Gyr, and 
used our input parameter a. When only cluster density is varied, 
rlmxb depends linearly on Y, however all the other parameters 
create significant scatter. We therefore urge not to overinterpret the 
dependence of the number of observed LMXBs on T, or the mean- 
ing of the power-law dependence of Y on density; the scatter in 
F in the observed GCs could easily be explained by the variance 
in other dynamical properties. If one assumes that in the observed 
sample of GCs the density increase is likely to be accompanied by 
increases of the cluster mass, and, accordingly, and a, then we 
expect from our results that the power law dependence of hlmxb 



ISivakoff et all d2007h . where metallicity effects are separated from 
dynamical effects. 

Even though a typical cluster contains 1 or 2 NS-MS LMXBs, 
there is only ~ 10 — 20 per cent chance that those LMXBs will 
be bright and persistent, and then only for metal-rich clusters. The 
probability to contain a bright persistent UCXB is no more than 
10 per cent and does not depend strongly on the metallicity. Over- 
all, there is at best a 30 per cent chance to contain a bright LMXB 
in a typical metal-rich cluster and a 10 per cent chance in a metal- 
poor cluster. This may explain the observed ratio between the bright 
LMXBs in metal-poor and metal-rich clusters. However, to confirm 
that, it is necessary to have observational statistics on qLMXBs in 
both metal-poor and metal-rich clusters - if the explanation above 
is correct, the number of qLMXBs will not show a dependence 
on the metallicity. Considering rich Galactic globular clusters, we 
find that our formation rate of UCXBs at 1 1 Gyr in Terzan 5 is 
about 10 per Gyr and in 47 Tuc, ~ 5 per Gyr (see Fig. @. Since 
tucxb ~ 10 s yr, this agrees with the presence of one bright (possi- 
bly ultracompact) LMXB in Terzan 5 and no bright UCXB detected 
in 47 Tuc. 

To compare the entire population of UCXBs in GCs with the 
results of our simulations, we integrated through the whole sam- 
ple of galact i c glob ular clusters, using data for velocities from 
iGnedin et alj J2002 PI and for other quantities from lHarrisI d 19961) . 
For this calculation, we used the formation rates shown in Table [6] 
and assumed that an UCXB will be bright for only ~ 10 s yr. As 
a result, the integral value of bright UCXBs in our Galaxy is 7.5, 
which is pretty close to the observed number of 4 to 9 UCXBs in 
GCs. The clusters that are most likely to contain UCXBs overlap 
very well with the clusters where those UCXBs are observed: NGC 
1851 and NGC 7078 each have a probability to contain an UCXB 
> 20 per cent. NGC 6712 is known to have been heavily tidally 
stripped, so its current cluster conditions do not re flect the con- 
ditions that produced its UCXB (see discussion in llvanova et alj 
l2005dh . Core-collapsed clusters are not properly modeled by our 
code, so our failure to predict UCXBs in NGC 6624 and NGC 6652 
(a likely UCXB) is not meaningful. Other clusters that have a > 
20 per cent chance to contain a UCXB include 47 Tuc, Terzan 5, 
NGC 6266, NGC 6388, NGC 6440, NGC 6441, NGC 6517, and 
NGC 6715 (M54). Unsurprisingly, this list consist of clusters with 
a high collision number. Those clusters have not been shown to 
contain bright UCXBs, but they generally contain numerous quies- 
cent LMXBs (some of which may be UCXBs) and MSPs (some of 
which may be products of UCXB evolution). Overall, it seems that 
physical collisions are able to produce UCXBs in numbers compa- 
rable to those observed in globular clusters in the Milky Way. 

It is more difficult to make a similar comparison for NS- 
MS LMXBs, as the lifetime at the bright stage is much less cer- 
tain. In the whole GC population we expect to have about 180 
qLMXBs formed through tidal captures and binary exchanges. We 
note that since our numbers for core-collapse clusters are lower 
limits, the real number of qLMXBs is expected to be higher. We 
identified those clusters that are predicted to contain 5 or more 
NS-MS LMXBs, most of which will be in quiescence. Those clus- 
ters are the same 10 as noted above for UCXBs: 47 Tuc, Terzan 
5, NGC 1851, NGC 6266, NGC 6388, NGC 6440, NGC 6441, 
NGC 6517, M 54, and NGC 7078. 5 of them, 47 Tuc, Terzan 5, 
NGC 6266, NG C 6388 and NGC 64 40 are known to contain 5 or 
more qLMXBs dHeinke et al"1l2003l . Cohn et al. 2008 in prep, ). 



with 7 is a bit less than one and close to the value of 0.8, as in 4 see on-line data at http://www.astro.lsa.umich.edu/5gnedin/gc/vesc.dat 
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NGC 6441 and NGC 6440 contain bright LMXB s with MS donors 
dVerbunt & Lewirj|2006l ; lAltamirano et alj|2007t) . We cannot make 
clear predictions for the core-collapsed, heavily reddened clusters 
Terzan 1, Liller 1, and Terzan 6 because we do not model dynami- 
cal evolution of clusters experiencing core collapse. We note that if 
AIC binaries contribute to LMXB formation, the numbers will go 
up, although not strongly as the characteristic MT time in this case 
is shorter than in the case of an average dynamically formed NS- 
MS binary. According to our models, several clusters have proba- 
bility 0.3 or more to contain a RG-NS LMXB, one of which - NGC 

7078 - indeed holds one. 

For extragalactic globular clusters, jBildsten & Delovell 20041) 
showed that the observed luminosity functions can be explained by 
a constant birthrate of UCXB binaries, 10 per Gyr per 200,000 Mq , 
for all GCs in a galaxy. We obtained almost constant birthrates for 
such binaries, but our rates are significantly lower (one per Gyr for 
200,000 Mq in a typical "dense" cluster). In fact, a 10 times higher 
UCXB formation rate will produce as many as 10 bright persistent 
LMXBs in Terzan 5. It is possible that the typical density of GCs in 
ellipticals is larger than that of GCs in our galaxy, which may ex- 
plain the variation in LMXB production. However, since we do not 
find any difference in UCXB production between metal-poor and 
metal-rich populations, our simulations do not support the idea that 
most X-ray binaries in extragalactic GCs are UCXBs. To explain 
most X-ray binaries in extragalactic GCs with UCXBs, we would 
need to assume that more NSs are retained (or formed) in metal-rich 
than metal -poor clusters (since the model of an irradiation-induced 
wind from lMaccarone et alj J2004I) is not applicable for degenerate 
donors). 



8 MILLISECOND PULSARS 

8.1 Populations of NSs that have gained mass 

In Table[8]we show the number of plausible pulsars that are formed 
in our simulations. In this table, we count as pulsars (i) all NSs that 
gained more than 0.01 Mq after their formation, through all the 
possible mechanisms of mass gain (mass transfer, common enve- 
lope accretion and during mergers or physical collisions) - includ- 
ing all recycled pulsars (see also fj3j ; (ii) all NSs that were formed 
through AIC/MIC within the last 10 8 years - "slow" high-B field 
pulsars (see also i]3,l| >. To investigate the role of mildly recycled 
post-AIC and post-collisional NSs, we include all of them in Ta- 
ble[8](c.f. EUand as well as all UCXBs that have evolved 

past 55 minutes (c.f. i )3.4t ; the reduction of the simulated popula- 
tions according to lifetime criteria will be done later in this Section. 
Three important facts can be seen from this table: (i) the numbers 
of ejected and retained pulsars are comparable; (ii) almost all NSs 
in binaries are pulsars; (iii) the numbers of produced pulsars are 
very large, more than several times the numbers of detected pulsars 
in specific comparable clusters (compare 320 MSPs in simulations 
vs. 22 detected MSPs for 47 Tuc, or 250 MSPs in simulations vs. 
33 detected MSPs in Terzan 5). 

The last point shows that a direct comparison of the numbers 
of NSs that gained mass with the observed numbers of radio MSPs 
fails, unless the detection probability is only about 10 per cent 
(however t his should be at leas t 33 per cent and may be even close 
to one, see lHeinke et al ]|2005bh . The alternative is that we produce 
and retain in our simulations rather more NSs than are retained in 
real GCs. In this case we cannot explain with our simulations the 
observed number of LMXBs. We conclude therefore that the num- 
ber of NSs that gain mass after their formation is correct, but, as 



was discussed in ij3] not all the mechanisms for mass gain lead to 
the formation of a radio MSP. 

Let us analyze the population of NSs which have gained mass 
in more detail. In Tables |9l and [Tol we show the results of simula- 
tions of 47 Tucanae and Terzan 5, distinguishing NSs with gained 
mass by their formation mechanism and by the process through 
which they acquired mass after the NS formation. Note that the 
plausible pulsars that were formed via AIC or MIC constitute more 
than half of all NSs that gained mass after their formation. As we 
discussed in §3.2\ such NSs may not reduce their magnetic field 
or may be failed to spun up and accordingly could have only very 
short life-times as a detectable MSP ("2a", "2b"). Indeed, very few 
of them were recycled recently, within last 10 8 years (see Tables|9] 
andllOt. From this moment, we remove from our further consider- 
ation of the population of "present" radio pulsars those AIC NSs 
that were recycled or formed more than 10 s yr ago (see also £13.11 

We produce MSPs i n AIC systems more efficiently than 
Kura nov & Postnovl J200^> . The most plausible explanation is the 
different adopted natal kicks distributions for ECS NSs, as they 
adopted a Maxwellian distribution with RMS 250 km/s (cf. 40 
km/s in our case). Nonetheless, we we find that AIC evolution 
does not lead to the formation of long-living radio MSP and there- 
fore can not explain most of the population of MSPs in GCs (c.f. 
Kura nov & Postnovl 12006]) ). 



8.2 Constraints from comparison with the observed 
populations of bMSPs 

In Figs.|6]and[7]we compare the simulated populations of NSs that 
have gained mass with the observations of 47 Tuc and Terzan 5. 
We use all simulated models for those clusters, and, as a result, the 
simulated populations plotted represent stellar populations about 2 
times larger than 47 Tuc and 5 times larger than Terzan 5. There are 
three panels for each cluster. Each panel shows the observed pop- 
ulation, compared to the simulated population of pulsars that were 
spun up via a different mechanism: (i) stable MT; (ii) CE or dynam- 
ical common envelope (DCE) during a physical collision with a gi- 
ant; and (iii) as a result of a merger. For most observed MSPs, the 
companion masses are calculated assuming a pulsar mass of 1.35 
solar masses and an inclination of 60 degree^], leading to some un- 
certainty in the direct comparison between observations and simu- 
lations. In the following we consider MSP populations, grouped by 
binary formation mechanism. 

8.2.1 Primordial bMSPs 

This population is present on the panels that show pulsars that gain 
mass via MT and via CE. The first population (MT) is character- 
ized by rather large periods, small companion masses, and low ec- 
centricities. Systems which may have been formed in this way exist 
in Terzan 5 (Ter 5 E), M53 (B 13 10+1 8), and M4 (B 16 20-26; note 
that a dynamical exchange is necessary to add the planet. lFord et al .1 
fcOOCh ). 

The second population, characterized by a companion mass 
around IMq (see the second panels for each cluster), is however 
undetected in both Ter 5 and 47 Tuc. No promising candidates are 
found among the entire GC bMSP population. Typical progenitors 

5 See http://www2.naic.edu~pfreire/GCpsr.html 
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Table 8. Plausible pulsars and their location at 1 1 Gyr. 
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The number of "plausible" pulsars that are retained in the halo or core of the cluster, as well as the number of ejected pulsars. As a pulsar we define here a NS 
that gained mass > 0.01 M . after its formation. "All psrs" - the number of all pulsars, "Bin psrs" - the number of binary pulsars and "NS-bin" - the number 
of all NSs in binaries, provided for comparison. For all GC models, except 47 Tuc and Terzan 5, the numbers are scaled per 200 000 Mq stellar population 
mass at the age of 1 1 Gyr; for 47 Tuc the numbers are given per its total mass taken as 10 6 Mq , for Terzan 5 - per 370 000 Mq . 



would be initially wide binaries containing two stars of rather sim- 
ilar masses from 6 to 10 Mq. Most of the primordial bMSPs are 
located in the halo. There are several possible explanations for these 
bMSPs to not be present in real clusters: 

• accretion during the CE either does not occur or does not lead 
to NS spinup and MSP formation; 

• as NSs in those bMSPs are mainly post-EIC (see Tables[9]and 
Hot - perhaps the NS natal kicks are higher, causing such systems 
to be destroyed rather than evolve via CE; 

• the CE prescription adopted in our code («ceA = 1) does not 
produce realistic results for these systems; 

• our mass segregation does not work properly for these systems 
- if the progenitors segregate into the cluster core faster, then they 
are destroyed and do not produce bMSPs. 

If indeed the kicks for EIC are higher than adopted in our code, 
then we will be able to retain very few NSs, as their main supply is 
from EIC. This will have a rather dramatic effect on the formation 
rate of LMXBs, decreasing it strongly. If the CE efficiency is the 
problem, then its reduction might help, as binaries will be formed 
with smaller binary separations, and will shrink their orbits faster 
and start MT. In our "CE-reduced" GC model (with CE efficiency 
is 10 times smaller) their number is a few times less, but still large 
compared to observations. Perhaps CE in these systems must be de- 
scribe d by a much smaller A parameter, e.g., |Pod siadlows ki et al.l 
d2003l) showed that A can be as low as 0.01. In this case all sys- 
tems evolved through a CE start MT or merge shortly after the CE. 
Alternatively, post-CE systems could be described by the angula r 
momentum loss prescription proposed in iNelemans & Toutl ( l2005h . 
Then post-CE systems will be much wider and thus more easily 
destroyed through dynamical encounters. Whatever the true rea- 
son, from now on we assume that in our simulations a CE event in 
primordial binaries with a NS formed via EIC does not lead to the 
formation of a MSP. 

Another population of systems which we do not show in Figs. 
6 and 7 are NSs produced in AIC, and then mildly recycled in the 
same binary system (see § 13,21 . An exampl e of a slow, likely post- 
AIC pulsar, is PSR B 1 7 1 8- 1 9 in NGC 6342 jLvne et alj 19931) . with 
a spin period of 1 s. Its location well outside the cluster half-mass 



radius suggests that it may be a halo object and primordial pulsar, 
perhaps produced through AIC. It is the only slow globular cluster 
pulsar with a companion in a circular orbit. Its companion's wind 
eclipses the pulsar, and the com panion may be filling its Roche lobe 
dJanssen & van Kerkwiikl2063) . which would support a recent AIC 
scenario. The cluster NGC 6342 is the least dense cluster inhabited 
by a slow pulsar, and PSR B 1718-19 has the shortest life-time (~ 
10 7 yrs) and highest B-field of any GC pulsar. Both characteristics 
suggest that this system may have been produced from a different 
mechanism than the other slow pulsars. The short life-time also 
implies a fairly high formation rate (see Tables 9 and 10). 

8.2.2 Dynamical common envelope events 

DCE leads at first to the formation of an eccentric binary MSP with 
orbital periods of 0.1-10 days, where a NS is spun up during the 
DCE. Later this system may start a new MT and become an UCXB. 
As a result, the population of DCE bMSPs can be found on two 
panels in Figs.|6]and|7]- among the NSs that gained mass via MT 
and via DCE. 

Binary MSP systems produced in NS-WD systems (from 
UCXBs) are produced in large numbers by our simulations, but no 
such MSPs have been observed in any globular cluster (see Figs. [6] 
and|7J. The possibilities for the fate of these systems are summa- 
rized in Q3.4\ We note here that failure of pulsar detection algo- 
rithms to detect their short periods seems an unlikely explanation, 
since the orbital periods of the majority of our simulated MSPs 
from UCXBs are similar to those of 47 Tuc R and Ter 5 A (see 
Figs. [6] and |7J. It seems that UCXBs do not produce observable 
MSPs in ultracompact binaries. In Table II II we retain statistics for 
these systems (as low-B, binary WDs with PjO. Id), although we do 
not know how they truly manifest themselves. 

The population of simulated bMSPs produced via DCE nicely 
explains observations of eccentric bMSPs, such as Ter 5 I, J, Q, 
X, Y, and Z, and 47 Tuc H. Our DCE simulations produce some 
binaries with low eccentricity, but cannot reach e < 10~ 3 . Thus, 
they cannot explain the low eccentricities of 47 Tuc E, Q, S, T, and 
U, as well as systems like NGC 6752 A. A stable mass-transfer 
mechanism must be invoked here (see discussion below in § 8.2.4). 
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Figure 6. Binary MSPs in 47 Tucanae. Observed and modeled populations: 
circles - observed, stars - formed via TC, squares - formed via DCE, tri- 
angles - formed via binary encounters and diamonds - "primordial" binary 
MSPs. Filled symbols represent binary MSPs with non-WD companions; 
in the case of observed systems filled symbols represent eclipsing bMSPs. 
Crosses mark eccentric bMSPs (e > 0.05), adjacent small rotated dia- 
monds mark bMSPs in the halo. The modeled population represents accu- 
mulated results from 5 models and corresponds to a stellar population twice 
as large as 47 Tucanae. 



Figure 7. Binary MSPs in Terzan 5. Observed (circles) and modeled pop- 
ulations, symbols as in Fig. [6] The modeled population represents accumu- 
lated results from 5 models and corresponds to a stellar population 5 times 
bigger than Terzan 5. 
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Table 9. Types of pulsars, at 1 1 Gyr, Terzan 5. 



Halo 

Total Binary 



Core 

Total Binary 



Table 10. Types of pulsars at 1 1 Gyr, 47 Tuc. 



Halo 

Total Binary 



Core 

Total Binary 



All "plausible" pulsars 



All 


56.1 ± 7.7 
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All "plausible" pulsars 
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The number of pulsars that are retained in the halo and the core of our 
Terzan 5 model. We indicate pulsars formed via CC, EIC, AIC or MIC, and 
that gained mass through different mechanisms (a mass transfer, a common 
envelope or a merger). 

8.2.3 Tidal captures. 

This population is present on all panels, as the MSP's spin-up gen- 
erally occurred before the MSP acquired its current partner ("5b" 
in Fig. [Q. The eclipsing bMSPs Terzan 5 P (P orb = 0.36 d, 
M c = 0.44M Q ) and 47 Tuc W (P orb = 0.13 d, M c = O.14M ) 
may have entered their current configurations through this route. 



8.2.4 Binary exchanges. 

The top panels of Figs. [6] and [7] include those bMSPs which were 
exchanged, and then underwent MT. While the companion mass re- 



Notations are as in Tabled but for our model of 47 Tuc. 

mains above ~ 0.2 A/0, the MT rate in these systems is high, and 
most of them will be seen as LMXBs or qLMXBs (see also £13 - 3b - 
However, the transition from LMXB to MSP is uncertain, and sev- 
eral bMSPs with nondegenerate companions are likely to have been 
formed through binary exchange or tidal capture (47 Tuc W and V; 
Ter 5 A, P, ad, and perhaps O; NGC 6397 A). Higher-mass bMSPs 
with nondegenerate companions may have their radio signals con- 
tinually obs£ured_by ionized material, leading to "hidden" bMSPs 
dTavanill 1 99 ll : lFreirell2005l) . 

Those bMSPs which underwent CV-like evolution after their 
exchanges can explain the very low-mass (brown dwarf) compan- 
ions to bMSPs (e.g. 47 Tuc I, J, O, P, and R; Ter 5 ae). CV-like 
evolution also can lead to the formation of binaries with longer or- 
bital periods (from one to ten days) and low-mass WD companions. 
We find that binaries with post-MT periods of 5-10 days are formed 
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relatively frequently, but all of them have collision times well be- 
low 1 Gyr. After the MT, they experienced a dynamical encounter 
that led either to an exchange of companions, or to an eccentric- 
ity increase. As a result, none of them are observed at 1 1 Gyr with 
low eccentricities. Potentially, we should be able to form binaries 
with a post-MT period of about a day. Such binaries have a col- 
lision time around 10 Gyr and thus should stay preserved in their 
post-MT shape for almost the entire cluster evolution. However, 
probably due to rather low number statistics, they are not formed 
in our simulations for 47 Tuc and Ter 5. Also, our post-exchange 
binaries are very eccentric. A better treatment of the MT in eccen- 
tric binaries, like discussed in ISepinskv et ail d2007h . may help to 
create low-eccentric binaries with low-mass WD companions. 

The bMSPS which were spun-up via mergers during a pre- 
vious dynamical interaction (bottom panels) generally have more 
massive companions and longer periods than the observed popula- 
tions. The few good candidates are Ter 5 Q, NGC 1851 A, M28 D, 
NGC 6441 A, and possibly M30 B. Dynamically formed binaries 
which undergo stable MT with a giant, producing non-eccentric 
long-period systems are rare in our simulations; these tend to come 
from the primordial population located outside the core. 

8.3 Refined populations of MSPs and observations 

From ij3. H and our comparisons with observations, we see that not 
all events where NSs gain mass lead to detectable radio MSPs. 
These events are likely channels for short-lived high- 1? MSP pro- 
duction: 

• recent NS formation via any ECS channel 

• AIC followed by continued accretion in the same system 

• mass gain via merger or physical collision 

Radio MSPs are not produced via these channels, or can not be 
detected yet: 

• accretion from a degenerate donor 

• mass transfer in an NS-MS systems with current MS mass > 
0.2Mq, as they have not yet turned on as radio pulsars 

• accretion during a CE event 

With these restrictions (see Fig. 1), we find that our estimates 
for the numbers of produced MSPs are 31 ± 6 in 47 Tuc and 26 ± 4 
MSPs in Terzan 5 (see Table ll It . These estimates are close to the 
observations for both 47 Tuc (22 detected MSPs) and Terzan 5 (33 
MSPs). 

Our simulations indicate that the fraction of single pulsars is 
higher in Terzan 5 than in 47 Tuc, as observed. The origin of iso- 
lated MSPs depends on the cluster dynamical properties. For ex- 
ample, in our standard model, about half of all isolated MSPs were 
produced by an evolutionary merger at the end of the MT from a 
MS companion ( t]3.4t . while the other half lost their companion 
due to a binary encounter (these MSPs were generally spun up by 
MT from a giant, leading to long periods). Most isolated pulsars in 
47 Tuc lost their companion due to a binary encounter; half of these 
involved the removal of a low-mass MS star at the end of its MT 
sequence. 

Spatially, MSPs are more likely to be located where they were 
formed - mostly in the core. Less than a third were produced in 
or ejected to the halo. The fraction of MSPs ejected into the halo 
increases as v ICC decreases. Roughly half of all o bserved GC pul- 
sars a re found outside their cluster core radius dCamilo & Rasiol 
|2005|) . but the radial distribution of most pulsars is exactly as ex- 
pected for a population that is produced in or around the core 



dGrindlav et al.ll2002l ; lHeinke et alj|2005bl) , with the ex ception of a 
few, l ikely ejected, pulsars in M15 and NGC 6752 (e.g. lColpi et al.l 
|2002|) . 

The creation rate of high-magnetic field (AIC) MSPs, at the 
age of 1 1 Gyr, is about 20-40 per Gyr in Terzan 5 and 47 Tuc. This 
leads to the probability of forming high-magnetic field MSPs (B < 
10 12 -gauss) with a period ^ 100 ms of less than a few per cent, 
while for B < 10 11 -gauss, it is about 50 per cent (s ee eq.[8). E.g., 
one p ulsar in Terzan 5 (J1748-2446J) has P=80 ms dRansom et al.l 
120051) . The absence of slow-period MSPs in 47 Tucanae suggests 
that the magnetic field in newly formed AIC NSs is most likely to 
be at the higher end, > 10 12 -gauss. We note that the results of our 
simulations are consistent with a higher probability of a high-B 
pulsar in Terzan 5 than in 47 Tucanae (see Table [TT), though the 
statistics are very poor. 

We compare the behavior of pulsar properties with increas- 
ing cluster density in our simulations, and in real clusters, in Table 
[T21 We divide the clusters containing MSPs with known proper- 
ties into categories roughly corresponding to the central densities 
used in the models: High-density (p c > 5.30, and core-collapsed 
with p c > 5.20): NGC 1851, NGC 6397, NGC 6522, NGC 6544, 
NGC 6624, M15, and M30. Terzan-5-like (5.05 < p c < 5.30): 
NGC 6266, Terzan 5, NGC 6440, and NGC 6441. 47Tuc-like 
(4.6 < p c < 5.05): 47 Tuc, NGC 6342, NGC 6626, and NGC 
6752. Standard (3.8 < p c < 4.6): M5, M4, and NGC 6760. 
Medium-density (2.8 < p c < 3.8): M53, M3, M13, NGC 6539, 
NGC 6749, and M71. No pulsars have been found to inhabit low- 
density (2.8 > p c ) clusters yet, nor do we produce any in our sim- 
ulations (Table [TT}. We exclude from consideration all MSPs with 
unknown binary properties, and those isolated pulsars in clusters 
which have been observed too rarely to measure binary properties 
(NGC 6624 D and E; NGC 6517 A and C) to reduce bias towards 
single MSPs. 

The column "All pulsars" indicates the total number of pul- 
sars used in that category; other columns refer to fractions of that 
number. The total numbers of pulsars cannot be usefully compared 
to the numbers of pulsars in the model, but the fractions of pulsar 
system types can be compared. Field systems in clude all field pul- 
sars w ith P < 100 ms and B < 5 x 10 10 G dManchester et al.l 
12005^ "BD" refers to brown dwarf systems; M CO mp < 0.07 M Q . 
"MS" identifies likely main-sequence companions, identified b y 
being eclipsing systems with M CO mp > 0.07 M© dFreirdl2~005l) . 
"He WDs" refer to all non-eclipsing systems with 0.07 < M c < 
0.5 Mq. "Heavy WDs" include all MSPs with M c > 0.5 M Q , 
unless the companion is known to be a NS. "> 10 d" refers to that 
subset of all MSPs with orbital periods longer than 10 days. "High- 
B" pulsars are taken to be MSPs with spin-period P > 100 ms. To 
improve the statistics for our high-B systems we count all high-B 
pulsars spun up within the past 1 Gyr, and divide this number by a 
factor of 10 to account for their shorter lifetimes. 

Several clear trends with the central density of globular clus- 
ters can be seen in the observations. Single MSPs rise smoothly 
from a rate of 0.25 of all MSPs in the field, to 0.7 of all MSPs in 
clusters. This is well-described by our model, suggesting that the 
formation of single MSPs may be reasonably well understood. The 
formation of MSPs with "brown dwarf" companions, with very low 
masses, exhibit a complicated density dependence, with few in the 
field, a quarter of all MSPs in standard-density clusters, and declin- 
ing fractions in higher density clusters. This indicates that their for- 
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Table 11. "Observable" pulsars at 1 1 Gyr. 



model low-B High-i? 

All psrs Single Binary All psrs 

Companion BD WD MS He WD heavy WD recently 

M c <0.5M Q >0.5M o recycled 

Period < O.ld < O.ld > lOd 0.1 + lOd 0.1 + lOd 

standard 5.2 ±1.7 2.1 ± 1.3 0.3 ± 0.3 1.0 ± 0.9 0.8 ± 0.5 0.5 ± 0.5 l-lij;! 0.5 ±0.5 0.4 ± 0.2 

metal-poor 3.9 ± 2.0 2.3 ± 1.2 0.2±|j | 0.5 ± 0.4 0.3±g | 0.2±° * 0.8 ± 0.6 0.2±° * 0.6 ± 0.2 

high-den 15.6 ± 2.7 9.8 ±3.1 0.5±g : | 3.5 ± 2.0 0.0 ± 0.0 1.4±J : | 3.0 ± 0.8 1-2±H 2.5 ± 0.5 

med-den 1.4 ±1.5 0.6±J g 0.0 ± 6.0 0.0 ± 0.0 0.3±g | 0.0 ± 0.0 0.2±0-:f 0.3±° : | 0.1 ±0.0 

low-den 0.0 ± 0.0 0.0 ± 0.0 0.0 ± 0.0 0.0 ± 0.0 0.0 ± 0.0 0.0 ± 0.0 0.0 ± 0.0 0.0 ± 0.0 0.0 ± 0.0 

low-cr 3.4 ±0.9 2.2 ±1.3 0.3±g g 0.3±g g 0.3±g g 0.7 ±0.7 0.2±g;| 0.0 ± 0.0 0.2 ±0.1 

long-t rh 1.6 ±0.8 0.8 ±0.6 0.2±g : ^ 0.3±g : | 0.3±g : | 0.0 ± 0.0 0.2±° : | 0.2±g-4 0.3 ±0.1 

BF05 3.0 ±1.5 0.9 ±1.0 0.1±[> : f 1.0 ± 0.4 0.7 ± 0.0 0.4 ± 0.4 0.3±g : ^ 0.3±g : | 0.3 ±0.1 

fast-MB 8.0 ±3.5 2.9 ± 1.6 1.1 ±0.9 1.5 ± 1.2 1.1 ±0.9 1.5 ±1.1 1.6 ±1.0 0.3±g : | 0.5 ±0.2 

CE-reduced 2.4 ±1.7 l-3±};§ 0.2±g;| 1.0 ± 0.4 0.3±°;| 0.2±g-| 0.5±g-g 0.2±° : 4 0.4 ± 0.2 

oldkicks 1.6 ±1.6 0.0 ± 0.0 0.3±g : | 0.0 ± 0.0 0.8 ± 0.8 0.5 ± 0.4 0.3±|] : ^ 0-2±oi 0.0 ±0.1 

47Tuc 25.7 ±3.5 9.4 ± 3.4 2.0 ± 1.4 6.1 ± 3.8 2.9 ± 2.7 1.6 ±1.7 7.8 ±2.7 2.9±| ; i 1.8 ±0.2 

Terzan5 30.8 ± 6.3 22.7 ±2.6 2.3 ± 1.1 9.1 ± 2.0 0.9±J;| 2.5 ±1.1 2.6 ± 2.3 0.8 ± 6.8 3.8 ± 0.3 

The numbers of pulsars that can be detected in a GC (including MSPs only if they gain mass via MT from a non-degenerate donor, or via DCE). Recent 
high-B pulsars are those which were formed or mildly recycled less than 10 8 years ago. Mass-transferring NS-MS binaries where a MS star is > 0.2Mq are 
excluded, as they are likely LMXBs/qLMXBs. "Companion" indicates the type of the companion; "M c " is the companion mass, and "Period" is the orbital 
period. For all GC models, except 47 Tuc and Terzan 5, the numbers are scaled per 200 000 Mq stellar population mass at the age of 1 1 Gyr; for 47 Tuc the 

numbers are given for 47 Tuc's mass, taken as 10 6 Mg; for Terzan 5, 370 000 Mq. 



Table 12. Detected pulsars and simulations. 



Cluster 


All pulsars 


Single 


BD 


He WD 


MS 


heavy WD 


> 10 d 


ecc 


high-B 


High-den model 
All high-den 


18.8 ±3.2 
19 


0.66 ±0.11 
0.71 


0.03 ± 0.03 
0.06 


0.16 ± 0.05 
0.0 


0.08±°;°1 
0.12 


0.10 ±0.08 
0.12 


0.01 ±0.00 
0.06 


0.19 ±0.10 
0.18 


0.14 ±0.04 
0.12 


Terzan 5 model 
All Ter5 


35.7 ±6.4 
49 


0.75 ± 0.06 
0.43 


0.05 ±0.03 
0.08 


0.08 ± 0.06 
0.31 


0.07 ±0.03 
0.12 


0.07 ± 0.03 
0.06 


0.03 ±0.03 
0.08 


0.12 ±0.07 
0.16 


0.11 ± 0.02 
0.04 


47Tuc model 
All 47Tuc 


28.3 ±3.1 
37 


0.38 ±0.06 
0.43 


0.07 ± 0.04 
0.19 


0.35 ± 0.09 
0.30 


0.07 ±0.06 
0.11 


0.14 ±0.12 



0.10 ±0.10 
0.03 


0.37 ±0.14 
0.08 


0.06 ± 0.01 
0.03 


Standard model 
all stand 


5.8 ± 1.7 
8 


0.44 ±0.18 
0.25 


0.04 ±0.04 
0.25 


0.27 ±0.17 
0.50 


0.12 ±0.10 
0.0 


0.15 ±0.10 



0.16 ±0.10 
0.25 


0.35 ±0.20 
0.25 


0.07 ±0.03 
0.0 


Med-den model 
All med-den 


1.5 ± 1.5 
11 


0.23±q;|| 
0.18 


n r)2+ - 04 
0.18 


n 27± - 39 
u - z,=c 0.27 

0.64 


n r)2+o.04 

u - uz ^0.02 




n 28+0-41 

U - ZO=C 0.28 




n nq+0.18 
0.18 


n 22+0- 39 
u - zz=c 0.22 

0.18 


0.07 ± 0.03 



Field 


77 


0.25 


0.04 


0.53 





0.13 


0.30 


0.01 





Fractions of GC pulsars by system type, simulations vs. observations for each central density. Observations summarized bv lCamilo & Rasiol HmO) . with 



additions from http://www.naic. edu^pfreire/GCpsr.html as of late 2007. See text (§8) for details. 



mation requires a high density environment, but also suggests that 
their binary evolution can be disrupted in high density conditions. 
Our simulations do not reproduce this trend very well; likely this is 
because our lack of a proper equation of state for brown dwarfs pre- 
vents us from accurately tracking the late evolution of NS-MS bi- 
naries and the eventual possible destruction of the companion. We 
may produce too many heavy WDs and too few He WDs, which 
both have only mild density dependences. The fraction of eccen- 
tric MSPs, perhaps surprisingly, shows little density dependence 
(in our simulations and in real systems). The production of MSPs 
with main-sequence companions and of high- 73 MSPs are sharply 
dependent on the c entral density of t he cluster, in agree ment with 
the suggestions by iFreird ( 12005) and iLvne et ail i ll 9961) that their 
production requires multiple binary interactions, and direct impact 



of NSs with other stars, respectively. Overall, our predicted num- 
bers are in quite good agreement with the observations. 



9 DOUBLE NEUTRON STARS IN CLUSTERS 

As predicted in we did not find primordial double neutron star 
formation to be efficient. In all our simulations - 70 cluster models 
or 2 x 10 7 Mq of stellar population evolved in different clusters 
- only 3 primordial DNSs were formed. In all three systems both 
NSs were formed via ECS, and none of the systems have merged 
within a Hubble time. 

Let us estimate the rate of dynamical DNS formation. The 
characteristic time for a single NS to have an encounter with a 
NS-binary can be estimated using eq. \10\ The binary compan- 
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ion in a NS-binary is usually less massive than the NS. If an ex- 
change occurs, the binary separation in the post-exchange binary a 
is larger than the binary separation ao in the pre-exchange binary, 
a ~ apmNs/m-2, w here m-2 is the mass of the replaced companion 
dHeggieetalJ 19961) . Binary eccentricities in hard binaries follow a 
therm al distribution, with an average eccentricity e ~ 2/3 l lHeggiej 
Il975t) . The maximum a that leads to a merger within 1 1 Gyr, if the 
post-exchange eccentricity is 2/3, is about 7.5Rq. The maximum 
mass for a NS-companion before the exchange is about the mass of 
a turn-off MS star - 0.9Mq . Then the maximum binary separation 
in a pre-exchange binary that leads to the formation of a merging 
DNS is a =4.8i? Q . 

From our simulations we find that the relative fraction of NS 
binaries in the core /NS-bm,core = AfNS-bin,corc/iVobjccts,core is 
0.175 per cent for the standard model, and 0.6 per cent for our 
Terzan 5 model. Then the characteristic time for a NS to have 
an encounter that may lead to the formation of a merging DNS is 
TDNS,m ^ 5 x 10 12 Gyr in the case of a standard cluster, in other 
words, about 500 NSs are needed to form one merging DNS within 
1 1 Gyr. This assumed that every encounter will lead to the exchange 
and that all NS-binaries have the separation of ~ 5Rq. In the case 
of Terzan 5, the chance of forming a merging DNS, per NS, is 25 
times bigger than in a standard cluster - the formation rate is highly 
dependent on the cluster density, as the power 2, and the number of 
retained NSs increases as the escape velocity increases. However, 
most NS-binaries have larger binary separations than ~ 5Rq, and 
even though encounters are more frequent in this case, they do not 
lead to the formation of a merging DNS directly. Subsequent en- 
counters of the wide DNS with other stars would be necessary to 
increase the DNS's eccentricity, and thus reduce its merging time. 

Only 14 DNSs are formed dynamically during 1 1 Gyr in all of 
our 70 cluster models. Since the formation rate of DNSs is highly 
density dependent, most of those DNSs were formed in the cluster 
models of Terzan 5. Only two of these 14 DNSs merged within 
1 1 Gyr, while another two DNSs merged within 14 Gyr, a Hubble 
time. 5 of the 14 DNSs stayed in the cluster for 1 1 Gyr, while the 
remaining 9 were destroyed within a Gyr after their formation. 

We note that the formation rate of DNSs is strongly con- 
nected to the formation rate of LMXBs, as the main building mate- 
rial of DNSs is NS-binaries that otherwise become LMXBs. As 
we show in ^7] we produce LMXBs in numbers comparable to 
those observed, but we underestimate the formation rates for core- 
collapsed clusters. In core-collapsed clusters it is especially clear 
that our simplified models do not properly treat the cluster dy- 
namics and therefore we cannot be very quantitative about the 
total number of DNSs that should be produced by all GCs, al- 
though only a small fraction of GCs are core-collapsed. We con- 
clude therefore that GCs are rather unimportant for the formation 
of merging DNSs and can at most provide a small (10 — 30%) 
contribution to the prod uction of short 7-ray bursts in old host 
galaxies as estimated bv lGrindlav et al .1 d2006l) . Such systems can 
only be produced in co re-collapsed clusters, l i ke Ml 5, where suc h 
a system is observed dAnderson et ai] ll990l : Ijacobv et aljExXih . 
Also, recent observations seem to indic ate that majority o f short 
GRBs are associated with young hosts dBerger et al.ll2007D . If in 
fact short hard GRBs o riginate from NS-NS mergers (although see 
iBelczvnski et alj|2007l) . then their rates and occurrence in young 
hosts is understood in terms of NS-NS formation without dynam- 
ical interactions, i.e. in field populations (Belczynski, Stanek & 
Fryer 2007). 



10 SUMMARY 

In our study of binaries with neutron stars in globular clusters we 
considered in detail the problem of the formation of LMXBs and 
MSPs. We predict from our simulations that most retained neu- 
tron stars in GCs must come from an electron capture supernova 
formation channel (c.f. the field, where most NSs are from core- 
collapse supernovae). A typical GC could contain at present (where 
the adopted cluster age is 1 1 Gyr) as many as ~ 220 NSs, of which 
half are located in the halo; a massive GC like 47 Tuc could have 
more than a thousand NSs. 

Analyzing encounters with NSs, we find that for NSs located 
in the core, about half formed a binary through an exchange en- 
counter and only a few per cent formed a binary through a physical 
collision with a giant or via tidal capture. The relative importance of 
tidal captures to physical collisions increases as the velocity disper- 
sion decreases. Although there are fewer binaries formed through 
DCE and TC than through binary encounters, the production of 
mass-transferring systems is roughly the same for all three chan- 
nels. We note as well that many physical collisions lead to mergers 
between NSs and other stars in the core. 

We derived a "collision number" from our globular cluster 
models which appears to scale linearly with the LMXB produc- 
tion frequency from our simulations, although only when the core 
density alone is varied. Variations of other GC dynamical properties 
with a fixed core density lead to a large scatter. This may explain the 
observed deviations from a linear dependence between the collision 
number and the number of LMXBs in non-core-collapsed clusters. 
We predict that the numbers of qLMXBs will have no dependence 
on metallicity, in contrast to the numbers of bright LMXBs. 

Our rates of LMXB formation predict 7.5 UCXBs in all 
galactic GCs, which is consistent with the observed number. For 
qLMXBs, we expect about 180 systems, whi ch agrees with the 
range (100-200) e s timated from observations dPoolev et alj|2003l 
iHeinke et al . 2003, 2005a). From our list of clusters which are ex- 
pected to have one or more LMXBs (bright or in quiescence) - 
47 Tuc, Terzan 5, NGC 1851, NGC 6266, NGC 6388, NGC 6440, 
NGC 6441, NGC 6517, M54, NGC 7078 - all clusters which have 
been studied with Chandra (all but NGC 6517) show at least one 
LMXB or qLMXB. Our results do not support the idea that the 
observed luminosity function in extragalactic globular clusters is 
principally due to UCXBs. 

The resulting retention fraction of ~ 5 — 7% (or ~ 220 NSs 
per 200,000 Mq), obtained in our simulations, seems to be firm, as 
the numerically derived formation rates of LMXBs, which strongly 
correlate with the retention fraction, are consistent with those ob- 
served. In our GC models, we achieved such a retention fraction 
by producing NSs via different channels of electron capture super- 
nova. If the natal kick veloci ty distribution will be revised in future 
from that used in this paper dHobbs et alj|2005h . then the necessity 
of ECS will decrease. If, however, it should be shown that AIC does 
not occur, then the retention fraction that is required to match the 
observations of LMXBs might have to be larger by a factor of 2. 

We find from our simulations that if all possible channels of 
NS formation and all possible mechanisms for their spin-up lead 
to MSP formation, then we overproduce MSPs. However, we still 
need these channels to produce observed LMXBs. We propose that 
high B-field MSPs (which are short-living) can be formed not only 
during core-collapse supernovae, but also due to physical collisions 
or accretion in a post-AIC system. We find that NSs which accrete 
and spin up during CE events overproduce bMSPs in the cluster ha- 
los from primordial binaries of intermediate masses. Such bMSPs 
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would be present in low-density clusters and have not yet been seen. 
In the case of the NS-WD LMXBs, we propose that the MT in such 
systems does not lead to the formation of radio bMSPs. The rates 
of UCXB formation (verified by observations of UCXB LMXBs) 
predict large numbers of ultracompact bMSPs in GCs which are 
not detected. 

Excluding the systems discussed above, as well as those which 
are still actively accreting their donor's material and are seen in- 
stead as LMXBs, we obtain lower limits of "detectable" bMSPs. 
The predicted numbers are in rather good agreement with the ob- 
servations - 26 ± 4 MSPs in 47 Tuc and 31 ± 6 in Terzan 5. The 
fraction of isolated pulsars is comparable to that observed and is 
larger in Terzan 5 than in 47 Tuc. 

Comparing the population census of our models with the ob- 
servations of all detected pulsars to date in GCs, we find good 
agreement for all types of pulsars - single, those with brown dwarf 
companions, with MS companions, with light or heavy WD com- 
panions, and with large periods. 

We do not find very efficient formation rates for double NSs 
- only a dozen were formed in 70 simulated models. These rates 
increase with the square of the GC core density. We conclude that 
DNS formation is most likely to occur in massive, very dense, and 
preferentially core-collapsed GCs, as suggested by the identifica- 
tion of only one double NS so far, in the core-collapsed GC M15. 

In conclusion we outline several important issues that must be 
addressed for further progress in studies of NSs in globular clusters: 

• Common envelope events will occur in primordial binaries of 
intermediate mass that produce NSs via EIC. What is the common 
envelope efficiency, and is the NS spun up by accreting the mate- 
rial? 

• What is the result of mass accretion onto a NS after it has ex- 
perienced a merger, either in a binary due to unstable mass transfer, 
or during a collision? 

• What is the final fate of a mass-transferring NS-WD binary? 

• How does the evolution of a mass-transferring NS-MS binary 
proceed when the companion's mass reaches < 0.04Mq? 

• What is the dependence of qLMXB numbers on the cluster 
metallicity? This holds the potential to determine the cause of the 
LMXB metallicity dependence (see § 7). 

In particular, if future observations show that the metallicity de- 
pendence of qLMXBs follows that for bright LMXBs, it will most 
likely indicate that these two types of clusters had different initial 
mass fun c tions (and thus star formation histories), as predicted by 
iGrindlavl i 19931) and expe cted from star formation theory in dense 
metal-rich environments dMurray||2007]) . If metal-poor GCs have 
similar numbers of MSPs as metal-rich GCs, but fewer qLMXBs 
(indicating a shorter LMXB lifetime), t hen the idea of irradiat ion- 
induced winds is likely to be correct dMaccarone et al. I l2004l) . If 
neither qLMXBs nor MSPs show a strong metallicity dependence, 
then th e idea of reduc ed magnetic braking in metal-poor GCs is 
correct (Ivanova 2006). 
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